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FOREWORD 


This  is  the  first  Technical  Report  of  the  Project 
on  Machine  Methods  of  Computation.  The  Project  Is  fi- 
nanced by  contract  N'jorlSo  of  the  Office  of  Naval  Re- 
search. It  is  an  outgrowth  of  the  activities  of  the 
Institute  Conunittee  on  Machine  Methods  of  Computation, 
established  In  November  1950* 

Purpose  of  the  Project  is  (1)  to  Integrate  the 
efforts  of  all  the  departments  and  groups  at  MIT  who  are 
working  with  modern  computing  machines  and  their  applica- 
tions and  (2)  to  train  men  In  the  use  of  these  machines 
for  computation  and  numerical  analysis. 
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INTRODUCTION 


This  study  deals  with  possible  lnnro\ naients  effected 
in  solutions  to  phj'sloc.l  nrobleTis  by  chenftes  in  the  coordinate 
system  of  the  probler.u  In  particalrr,  perturbation  solutions 
to  the  problems  are  Investigated  with  a view  to  Improving  the 
usual  solution  by  allowing  a pei-turbatlon  of  the  Independent 
variables  s.s  uell  as  of  the  dependent  variables.  The  technique 
has  proved  useful  foi’  wave  propagation  problems  of  a hyperbolic 
nature  described  in  Chapter  2,  but  has  been  less  successful 
for  the  elliptic  problem  of  the  flow  past  a thin  airfoil  dis- 
cussed in  Chapter  3. 

The  idea  of  perturbing  the  Independent  variable  of 
a problem  was  used  some  time  ago  by  Poincare  in  finding  the 
limit  cycle  of  a nonlinear  oscillator  , and  more  recently 
(1949)  M,  J.  Llghthlll  considered  the  method  at  some  length 
In  a most  Interesting  paper  (12)  , Llghthlll  considers  vari- 
ous types  of  ordinary  differential  equations  whose  usual  per- 
tiirbatlon  solutions  are  unsatisfactory,  and  also  In  this  psper 
he  deals  briefly  with  nonlinear  partial  differential  equations. 
In  a later  paper  C^l)  he  applies  the  technique  of  coordinate 
perturbation  to  correcting  the  usual  solution  to  the  problem 
of  Incompressible  flow  past  a thin  airfoil.  In  the  usual 
solution  which  Is  based  on  the  small  slope  of  the  airfoil 
profile  a satisfactory  answer  is  not  obtained  near  the  lead- 
ing edge  of  a blunt  nosed  airfoil,  but  Llghthlll  found  that 
a small  constant  shift  of  the  coordinate  system  corrected  the 
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firat  order  solution,  luil-lng  It  valid  uniformly,  and  he 
Implied  th&t  Bucoeeelve  oorreotions  could  be  made  in  the 
coordinates  as  the  higher  order  perturbation  solutions 
were  found. 

This  problem  vas  the  initial  impetus  of  our  in- 
vestigation, and  it  was  planned  to  study  the  problem  at 
greater  length  and  in  particular  vith  reference  to  the  case 
of  oompreeslble  flow.  However,  as  the  study  progressed  it 
became  clear  that  the  perturbation  series  developed  for  the 
Independent  variables  were  not  satisfactory.  Actually  the 
Bucceseivo  perturbation  f^onotlons  in  these  series.  In  order 
to  correct  the  increasing  order  of  singularities  In  the 
velocity  functions  near  the  leading  edge,  have  to  carry  the 
singularities  themselves,  and  in  the  singular  region  each 
term  of  the  perturbation  series  becomes  of  the  same  instead 
of  decreasing  order.  The  results  and  the  failure  of  the 
series  near  the  leading  edge  are  discussed  in  more  detail 
in  Chapter  3» 

Nevertheless  the  hope  remained  that  at  least  some 
sort  of  one-stage  correction  could  be  found  for  the  case  of 
compressible  flow  past  a thin  airfoil,  since  the  small  con- 
stant coordinate  shift  does  effectively  correct  the  Inoom- 
presslble  case.  Several  lines  of  investigation  were  attempted 
and  various  types  of  coordinate  correction  were  considered, 
but  iinfortunately  none  of  them  proved  successful.  For  example 
a constant  coordinate  shift  Is  no  longer  correct  since  it  dues 
not  9atl?;fjr  the  more  complicated  equations  for  the  compressible 
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case.  Then  as  far  as  other  types  of  functions  are  concerned 
one  Important  aspect  arising  from  the  Investigation  should  be 


rcma:(*ked.  Upon  attempting  series  exi’'an8lons  baaed  on  a per- 
turbation parameter,  £.  , In  the  form 

X = x^°^(X,Y)-h  + L^x^^^vX.Y)  ^ (1.1) 

where  X and  Y are  new  coordinates,  It  was  found  that  functions 

/ 4 \ 

of  the  type,  , appearing  In  (1.1)  ooald  not  be  obtained, 

but  rather  chat  due  both  to  boundary  conditions  and  to  changes 
of  order  In  the  derivatives  of  the  perturbation  functions  that 
the  functions  were  of  the  form 

x^^^(X,Y,6)  (1.2) 

In  other  words  the  functions,  x^^\  did  not  fit  the  usual  type 
of  perturbation  solution  (1.1). 

Actually  such  functions  might  be  acceptable  In  a new 
frame  of  perturbation  solutions  of  the  t3rpe 

X s x^°^(X,Y,t)+ e.x^^^(X,Y,6)  + (X,Y,£)  + ...  (1.3) 

where  the  £ necessary  to  fulfill  the  requlrsaents  of  the 
problem  Is  allowed  to  appear  In  the  perturbation  functions. 

This  line  of  Investigation  has  not  been  pursued,  but  might 
lead  to  new  methods  of  problem  solution  and  Is  perhaps  worthy 
of  some  study. 

Next,  since  It  had  not  been  possible  to  find  real 
solution  Improvement  for  the  elliptic  case,  the  Investigation 
turned  to  problems  of  a hyperbolic  character.  In  particular 
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the  nonlinear  hyperbolic  partial  differential  equations  govern- 
ing wave  propagation  dependent  on  one  space  variable  and  one 
time  variable  were  oonaldered.  Problems  of  this  sort  for  the 
oases  of  spherical  and  oyllndrloal  waves  have  been  discussed 
by  Whltham  (23) , (24) , in  papers  where  a correction  In  one 
set  of  ohauraoterlstlos  Is  advocated.  That  Is  In  lieu  of 
using  the  usual  functions  of  a linearized  characteristic  vari- 
able In  the  form 


t{x  - c^t)  , 


c^  s speed  of  sound  at 
Infinity 


(1.4) 


where  the  characteristics  ©j?e  straight  and  parallel  in  the(x,  t) 
plane,  Whltham  suggests  a correction  In  the  form  of  a new 
variable,  z,  which  Is  to  be  constant  along  the  true  character- 
istics of  the  problem  now  to  be  solved  by  functions  of  the 
type  f(z).  He  h-as  obtained  interesting  results  from  this  ap- 
proach, but  a drawback  of  the  method  is  Its  failure  to  deal 
with  both  sets  of  characteristics  at  the  saine  time. 

The  question  arose  as  to  whether  It  might  not  be  pos- 
sible to  devise  a method  correcting  the  solution  along  both 
characteristic  directions,  and  thus  sllow  for  waves  propagating 
In  either  direction.  It  seemed  natural  to  consider  solutions 
exoanded  In  terms  of  characteristic  variables  both  for  the 
dependent  and  Independent  variables . since  In  this  way  both 
the  physical  quantities  and  the  correct  characteristic  curves 
might  be  found  In  the  physical  (x,t) -plane  by  a mapping  from 
the  characteristic  plane.  Attention  was  restricted  iralnly  to 
the  case  of  plane  wave  propagation  where  the  equations  are 


simpler  than  for  higher  dimensional  csees,  and  a solution  was 
attempted  in  terms  of  charrcteristlc  yarieblee  o<  and  ^ for 
velocity,  u,  and  spetul  of  sound,  c,  and  for  the  originally 
Independent  space  and  time  variables  x and  t.  For  a per- 
turbation parameter,  £.  , the  solutions  were  assumed  in  the 
form 

u = u^°^(o(,p)4-  £ u^^^  (C3<,p)  -f-  4-... 

c « (cK.,  p)4-  £c^^^(oc,p)  -h  4-  ... 

X = x^°^(oi,p)4-  4-  4-  ... 

t = t^®^(o^,p)4-  £t^^^  (o<  ,p)  4-  {«v  ,p)  4-  ... 


(1.5) 


Such  a typo  of  solution  proved  to  be  most  satisfactory. 

In  the  case  of  a plane  wave  the  series  for  u and  c terminate, 

and  a general  expression  can  be  found  for  the  nth  perturbation 

function  in  the  series  for  x and  for  t.  Furthermore,  and  of 

some  Importance,  one  finds  that  the  convergence  of  these  latter 

two  earles  can  be  demonstrated  for  sufficiently  small  Initial 

disturbances  so  that  the  correctness  of  the  solution  Is  assured. 

More  precisely,  the  convergence  is  dependent  on  the  slae  of  the 

Initial  disturbance  u'  and  o'  to  the  base  quantities  u^  and 

o o 


only  to  the  extent  of  requiring  that 


u'  ✓ 


2 

TT 


o'  ^ 1 


(1.6) 


where  X Is  the  adiabatic  exponent  so  that  for  example  If 

y = 1-^ 


ul  < 0.833 


c ' / 0,166 


I 
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In  terme  of  density,  ^ , this  becomes 

^ < 1.161 

which  le  of  the  order  of  the  requirement  on  velocity  dle- 
burbance.  In  other  words  the  series  is  a convergent  one  even 
for  an  Initial  density  ratio  as  high  as  2.16, 

The  convergence  of  the  series  in  the  ordinary  regions 
of  the  flow  is  useful,  but  the  outstanding  fact  is  that  conver- 
gence persists  into  the  region  where  the  mapping  from  the  charac- 
teristic plane  onto  the  physical  plane  becomes  multiple-valued 
and.  where  a shock  would  develop.  In  such  a region  the  charac- 
teristics in  the  physical  plane  form  an  envelope,  and  one  can 
show  from  the  solutions  (1,5)  that  the  envelope  develops  for  t 
of  the  order  of  1/6  • Then  by  digressing  to  state  that  the 
usual  type  of  perturbation  solution  takes  the  form  (see  Chapter  2) 

u - ^®\x,t)-l-  £f  (x,t) E^tf  (x,t) -f  £^t  f^^^(x,t)-t-  ... 
or  (1.7) 

u = Mx,t)+  £ j^f^^^(x,t)+  fctf^^^(x,t) -f  (£.t)^f^^^(x,t)-r  ...| 

the  great  improvement  introduced  by  solution  (1,5)  becomes  olaar. 
Equation  (1.5)  converges  in  a perturbation  sense  of  decreasing 
order  of  terms  even  as  t becomes  0(1/6.)*  out  the  terms  In  (1.7) 
become  all  of  the  same  order  for  this  size  t and  no  longer  repre- 
sent a satisfactory  perturbation  solution.  The  solution  (1.7) 
falls  in  t,h<>  some  v^ay  as  that  for  thf  elliptic  case  did,  but 
(1.5)  r-5mai:if  actor?. . 
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In  extending  (1,5)  into  the  region  where  the  solution 
hecomefl  multlple-vnluod,  a shock  wave  is  introduced  to  allow  a 
Jump  In  the  phyaioal  quantities.  Its  location  can  be  satis- 
factorily demonstrated,  and  one  can  show  thet  the  shock  speed 
is  the  rverege  of  the  slopes  of  the  characterlatlcs  on  either 
side.  The  distance,  x,  which  the  solution  may  be  continued 
beyond  the  point  of  shook  development  depends  on  the  degree 
of  accuracy  desired.  It  Is  shown  that  an  error  of  the  order 
of  Is  Introduced  into  the  solution  at  a distance  x along 

the  shock. 

The  theorj  developed  In  Chapter  2 is  Illustrated 
there  by  a particular  example  of  an  Initial  value  problem  In- 
volving a periodic  distribution  of  density  In  a fluid  Initially 
at  rest.  The  solution  to  the  problem  Is  found  both  for  the 
usual  type  of  perturbation  theory  and  for  the  Improved  type, 
and  the  advantages  of  the  second  type  are  demonstrated. 

One  would  hope  that  the  same  advantages  arising  from 
coordinate  perturbation  In  the  plane  wave  case  would  carry  over 
to  the  higher  dimensional  cases  of  cylindrical  and  spherical 
flow.  Of  course  one  can  carry  out  the  same  steps  of  the  process, 
but  for  these  oases  where  the  equations  are  more  complicated  the 
series  for  u and  c no  longer  terminate,  and  as  yet  no  oonveigence 
proof  for  the  series  has  been  obtained.  There  Is  thus  no  guarantee 
of  the  solutions,  but  as  Is  after  all  usually  the  case  with  per- 
turbation solutions  one  can  assume  validity  for  •small  enough* 
epsilon.  It  seems  appropriate  to  extend  the  coordinate  per- 
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turbation  technique  to  theue  problens,  Justifying  the  procedure 
by  the  improvement  the  technique  brought  to  the  plane  wave 
oaee,  and  hoping  tbat  In  time  a convergence  proof  will  be 
developed  outlining  the  regions  of  validity  of  a solution. 

In  summary  then  coordinate  perturbation  Is  a very 
useful  technique  for  Improving  the  perturbation  solutions  to 
problems  of  a hyperbolic  nature.  The  oharacterlstlc  variables 
are  the  "natural"  variables  of  the  problem  and  It  now  seems 
clear  that  all  the  quantities  of  a problem  Including  the  phyei- 
cal  coordinates  should  be  expanded  In  terms  of  the  character- 
istic variables  If  the  best  type  of  perturbation  scheme  Is  to 
be  used.  For  the  elliptic  case  on  the  other  hand  no  sot  of 
natural  variables  Is  available.  Problem  solution  In  terms  of 
the  imaginary  characteristics  Is  of  no  advantage,  and  so  when 
a solution  becomes  unsatlsf aotory  and  slngalarltlee  arise 
ore's  only  recourse  Is  to  modify  their  effect  In  some  way.  In 
the  case  of  the  airfoil  the  singularity  could  be  partially 
hidden  Inside  the  profile  but  Its  presence  was  still  felt  In 
the  region  of  the  leading  edge  end  the  singular  behavior  of 
the  solution  in  the  region  persisted. 

It  la  concluded  that  coordinate  perturbation  as  dis- 
cussed in  this  thesis  does  not  yet  offer  a real  Improvement  to 
elliptic  problems,  but  rather  finds  great  usefulness  In  Its 
application  to  hyperbolic  problems  where  the  technique  allovrs 
the  solution  to  hold  even  somewhat  after  the  development  of  a 


shock. 
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THE  TOERBOLIC  CASE 

A.  Plfine  Vave  Propag-atloa 
Introduction: 

The  solution  to  the  problen  of  one-diinenelonnl 
Isentroplc  vave  propagation  may  alv’ays  be  given  In  the  form 
of  an  Integral  Involving  the  Rlemann  function  and  the  Initial 
conditions  and  boundary  conditions  of  the  problem  (16]) . 

However,  except  In  certain  special  cases,  such  a representa- 
tion of  the  solution  Is  not  very  tractable,  and  some  other 
form  would  be  more  useful. 

One  of  the  alternative  forms  of  solution  often  ap- 
propriate, and  In  fact  one  of  the  favorite  approaches  of  applied 
mathematicians.  Is  that  of  a perturbation  solution  where  the 
problem  is  linearized  by  expanding  the  solution  in  terms  of 
a small  parameter  of  the  problem.  The  resultant  solution  Is 
satisfactory  as  long  as  the  series  beha.ves  properly,  but  there 
may  be  regions  where  the  solution  falls.  In  peo'tlcular,  as 
will  be  shown  later,  for  the  problem  of  plane  wave  propagation 
the  series  will  diverge  In  regions  vrhere  a shock  starts  to  form 
and  In  these  regions  the  series  solution  Is  no  longer  satisfactory. 

In  the  attempt  to  overcome  such  drawbacks  of  the  per- 
turbation t3Hpe  of  solution,  and  in  line  with  the  general  alms 
of  this  Investigation,  a solution  to  the  plane  wave  problem 
has  been  found  In  terms  of  a perturbation  not  only  of  the 
dependent  quantities,  velocity  and  density,  but  also  of  the 
space  and  time  variables,  x and  t.  The  nevi  Independent  variables 
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used  are  the  charaoterlatio  parameters,  and  the  four  perturba- 
tion aeries  for  x,  t,  velocity,  and  density  are  expanded  In 
terms  of  functlor®  of  the  characteristics , The  solution  seems 
most  promising  In  that  the  convergence  of  the  series  can  be 
proved  for  an  appropriately  chosen  perturbation  pso*ameter,  and 
also  In  the  fact  that  the  region  of  ahoolc  development  can  be 
studied.  That  Is,  even  though  the  mapping  of  the  characteristic 
plane  onto  the  physical  plane  becomes  multiple-valued  indicating 
the  appearance  of  a shock,  the  solution  series  does  not  break 
down,  and  may  be  used  to  describe  the  phenomenon  even  somewhat 
after  the  formation  of  a weak  shook. 

In  the  first  section  of  this  ppjt  the  four  perturbation 
series  are  found  for  the  initial  value  problem  with  general 
Initial  conditions,  end  a proof  of  the  convergence  of  the  series 
Is  given.  An  example  of  a simple  physical  problem  is  then  given 
to  Illustrate  such  aspects  is  the  Inadequacies  of  the  usual  type 
of  perturbation  solution  and  the  Improvement  effected  by  the  Im- 
proved perturbation  solution.  A further  simplification  of  the 
example  problem  is  then  made  so  that  the  solution  may  be  found 
easily  from  the  Rlemann  function,  expanded  In  terms  of  the  small 
perturbation  parameter,  g,  , and  compared  with  the  perturbation 
solution.  The  solutions  are  found  to  agree  and  in  fact  the  order 
cf  £.  required  to  allow  the  expansion  of  this  latter  solution  Is 
that  expected  from  the  general  theory. 

The  second  section  Is  concerned  with  the  Incidence  of 
a shock  wave.  The  geometry  of  the  mapnlng  of  the  characteristic 
plane  onto  the  physical  plane  Is  discussed,  and  the  ability  of 
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the  solution  to  penetrate  ellphtly  Into  the  region  where  the 
mapping  starts  to  "fold  over"  Is  demonstrated.  In  the  dis- 
cussion entropy  changes  are  neglected  so  that  acctiracy  only 
through  second  ci*ler  in  sliock  etrength  Is  maintained.  The 
same  physical  example  used  cbove  Is  described  near  the  regions 
of  shock  development. 

Section  1;  Perturbation  solution  and  proof  of  convergence. 

When  the  equations  governing  the  propagation  of  a 
plane  vave  are  cast  into  the  form  of  differential  equations 
along  the  characteristic  directions  of  the  problem,  they 
become,  for  the  case  of  a polytroplc  gas, 


X,,,  ■ (u  + o)t*< 
a«  (u  - o)t^ 


(2.1) 


where  )(  is  the  adiabatic  exponent,  o<  and  ^ , the  character- 
istic variables,  and  u and  o the  local  velocity  and  speed  of 
sound  respectively  (^4^. 

Since  a perturbation  form  of  the  solution  to  the 
initial  value  problem  Is  to  be  obtained,  the  initial  conditions 
must  be  given  ae  disturbances  superlmoosed  on  some  prior  uniform 
state.  The  choice  of  values  to  be  assigned  to  q<  and  ^ on  the 
Initial  line  Is  governed  only  by  the  requirement  that  a correct 
parametric  (for  example  single-valued)  representation  of  the 
Initial  state  Is  achieved  since  the  equations  (2.1)  are  homo- 
geneoub  in  or  ^ . For  simplicity  the  initial  line  t « 0, 

in  the  characteristic  plane  will  be  chosen  as  x » et  « ^ . 
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Thon  for  a perturbation  paraaeter  £ , euid  Initial  perturbations 
and  o^tix)  to  o and  u,  the  Ixiltlal  oondltlone  beccsas 

^ X » ei  « ^ 

t » 0 

“*?  ■(  c = 0^(1  + tg(«))  = e„(X+  Sg(p)) 

(^U  « eCQf(«)  a £0^f(^), 

where  the  unperturbed  gae  was  aseuaed  at  rest  with  uniform 
density. 

A solution  to  equations  (2.1)  under  the  conditions 
of  (2,2)  Is  to  be  attempted  in  the  form  of  the  perturbation 
series, 

X - x^°^+  £x^^^(o(,^)  4-  e^x^^^(«,«)4-  ... 

t s t^°^+  ct^^^(c^,p)  -h  (o<,^)4-  ... 

(2.3) 

u » u^®^-!-  £.u^^^(ot,^)  £^u^^^  («,^)^  ... 

c = £c^^^(o<,^)  4-  t^c^^^(rt,p)4-  ... 


This  form  of  solution  will  linearize  the  problem  and  give  seta 
of  equations  from  which  the  higher  order  perturbation  functions 
may  be  successively  determined. 

The  terms  free  from  £ In  the  equations  and  Initial 


conditions  yield  the  Initial  approxlmatlcn 
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.(0)  o(-i-g 

= 2^ 

(0)  a - g 


= 0 


- c 

c - Cq. 

For  the  functions  u and  c,  the  boundary  conditions 
together  with  the  two  equations  of  (2.1)  homogeneous  In  u and 
c,  show  that  In  (2,3)  ^ne  will  obtain 

u^^^  = 0 and  =0  for  k > 1,  (i 

whereas  for  u^^^  and  one  finds 

li 4-2-—  - r(B) 

2 X -1  V 

(: 

,,(1)  ^(1) 

u C 

J. 

2 X -1 


where  the  Initial  conditions  require 


*•<?>  = 

\ 2 U "J- 

-“■••(‘f -IS 


Prom  these  expressions  and  from  equations  (2.1),  any 
(ki  (ki 

perturbation  function  x'  ' or  t'  ' Is  found  from, 


J 
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X = 


ov  . 

U +-C 


(K-0 


(2.7) 


In  equations  (2,7)  using  equations  (2,6)  the  Integrands 


contain. 


2c^ 


s Ar(p)  - Bs(pt) 


(2.8) 


where 


t 4-  1 


A = 


• = Br(^)  - As(o<) 

3 - ^ 


a ^ 


0 o 


(k ) (k^ 

Now  since  the  functions  x'  • and  t'  • ai»e  analogous,  only  the 
(k) 

functions  t'  ' will  be  treated,  and  It  will  be  shown  that  these 
functions  may  be  expressed  In  a certain  general  form.  Then  from 
this  expression,  the  convergence  of  the  perturbation  series  for 
t may  be  proved. 

The  form  assumed  by  the  functions  will  be  showi  to  be. 


c.  t'”*'’  = ^ r'""  r'(f)  S'(.)  j r '(5)  ^ 


(2.9) 


where 


15. 


^ ?,v 


r R 


(n+D^’^ 


= constants  such  th^it: 


b!; 


s(h') 


(2.10) 


^ - kronecker  delta 
>o 

(all  indices  in  (2.10)  non-negative) 


Proof  of  the  validity  of  (2.5)  Is  carried  out  by 
mathematical  Induction  from  equations  (3*7)  (2.8). 

At  the  first  stage  for  = 0 

evidently  (2.9)  holds  since 


(0)  _ & 


From  (2.9) : 


•T>l 


V 


r 


-j(2.1l) 


so  that  equation  (2.7)  for  t^ becomes 


and  equation  (2,12)  may  "be  Integrated  using  (2.13)  tc  give 


(3.1^) 
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This  expression  la  of  the  form  (2.9)  and  upon 
proper  treatment  of  Indices  yields  the  rnnurrence  forrauli, 
(2.10).  From  (2.10)  of  course  one  finds  that  the  rule  for 
symmetry  In  the  Indices  (by  pairs  -U,m)  and  (ju,y))  holds 
at  the  (n+1)—  stf  ro  If  It  holds  at  the  n— . 


The  next  step  In  the  lnvestlf;itlon  is  the  proof  of 
the  convergence  of  the  series  for  t, 


... 


from  the  expression  (2.5)  for  the  functions  t 


(k) 


(2.15) 


It  will  be  shown  that  If  the  functions  r and  s are 


bounded  for  the  entire  range  of  their  argument  by,  say,  a 
constant,  M, 


»rl  ^ M 

|a|  4 M, 

(which  means  of  course  from  (2.6b)  that  the  original  per- 
turbation functions  f ajnd  g must  be  bounded)  , then  the  fol- 
lowing geometric  aeries  may  be  shown  to  dominate  (2.15), 


(2.16) 


(2.17) 


Tills  series  converges  If  g,  Is  chosen  such  that 


£<7^  . 

4am 


(2.18) 


For-  proof,  the  terras  In  the  series  for  o t are 

o 


evaluated  from  (2.9): 


18 


t 


L Nr («-p)  1 


or 


(2.19) 


n*»' 

C.1  i 


4 1-“  M-  (.-p)  21 1 r“  ;t,  1 


n+t 


The  problem  then  resolves  Itself  Into  estimating 
the  exprsEsion 


f,v 


(2.20) 


The  estinjntlor.  Is  nrde  by  siur  nlng  the  recurrence 
formula  (2.10)  In  terms  of  absolute  values.  J\irthermore  of 
course  only  three  of  the  four  Indices  ^ P • ^ 

/ I 

Independent  so  one  may  set 

n+i- . (2.21) 

Finally,  to  make  use  of  the  delta  functions,  the  sum  is  broken 


up  into  the  form 

n-t-1  n-»-l-in  n*-l  m n-*-l-m 

21  21  z = z z:  z: 

maO  p+v  am  jXaO  msO  y«0  *0 

where  ^ z n - V » 

which  represents  correctly  the  requirement 

0-^  ^ n+1,  where  no  arrangement  is 


(2.22) 


fr 


repeated. * 


Then  from  (?,.10) 


nM-m , . . 


n*l  m I . .1 

- Z ZZ|  Rir;'  1^''^- a)  + K'zz  j 

m-a  ^‘■o\  ’ ’ / 


■'I 


("1 


■yy  T h-l- 

lL{  L I V'  ■'  „ 7 

“ m, o ( 

* |»o 


(M 


Dq4i  ' A ” 

b 


\ 


I 

/ 


(2.:>3) 


In  thl8  expression,  from.  (?..8) 


A = 


V4-  T 

w * A. 

c 

1^ 

U’Q 


so  that  for  If  > 1,  A > B and 


- A ^ A . 


(2.24) 


n+i 


* To  prove  that  .„  , , 

III  r“ 

does  Include  every  type  of  R,  It  suffices  to  show  that  s partlculoj 
R (l,J,k,)!  = definite  values  such  that  1+-J  k » n+1) 

can  appeeir  In  one  and  only  one  way«  In  other  v'ords 


m can  only  equal  k 

* * * i 

and  « « " J 


vim  as  required, 
^in+l-m  as  required. 
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Then,  by  replacing 


by  1 In  (2.23),  and  by  showing 


q 1 

that  In  each  of  the  form  sums,  expressions  of  the  type, 


z 

(the  proof  In  each  cese  Is  the  sane  as  that  at  the  bottom  of 
page  19).  one  acbdevea  an  evaluation  of  (2.20) 


(n4:>  > 


^ I R' ' ;r. 

(n) 


A-^  >A 

I 


4A2_ir 


A-v^-^p  + V 


(2.25) 


th 


So  that  If  T * Is  an  upper  bound  at  the  n stage, 


y I ^ J 4 4AT 

£_  \ n K y ’ ■< 

A4^*P4>.  0 + 1 

Finally  then,  from  (2.19), 


(2.26) 


SO  that  (2.17)  does  represent  a dominant  series  for  (2.13), 
requiring  a restriction  of  £,  of  the  form  (2.18)  for  con- 
vergence, and  the  proof  Is  completed. 

Illustrative  example 

The  points  discussed  so  far  In  this  section  will  be 
illustrated  with  reference  to  a particular  simple  example  where 
the  function  f(x)  of  equation  (2.2)  Is  taken  to  be  zero,  and 
the  function  g(x)  to  be  cosx. 


Thus  at  t 


0 


’ u = 0 

e = Cjjd  4-  £ coax) 


(2.27) 


*N 


1 


I 


I 


i 

1 


The  general  behavior  of  an  Initial  periodic  dis- 
turbance of  this  port  such  as,  for  example,  the  development 
of  the  higher  order  hermonics  and  later  the  tendency  of  shockr 
to  appear  Is  of  course  of  some  Interest  In  Itself. 

1,  Usual  perturbation  aolxtllon 

In  the  Introduction  It  was  mentioned  that  the  usual 
type  of  perturbation  solution  In  which  only  the  dependent 
quantities  u and  c are  expressed  as  perturbation  series  often 
Is  unsatisfactory.  The  particular  way  In  which  the  solution 
Is  Inadecuate  for  this  exeunple  problem  will  be  demonstrated 
below. 

The  equations  for  plane  wave  propagation  referred 
to  X and  t as  Independent  variables  are 

y - 1 

c^  + uc^  H cu_  B 0 

(2.28) 

cc^  + — + uu^)  = 0. 


In  the  usual  pertxirbatlon  scheme  a solution  to  (2.28) 
subject  to  the  initial  conditions  of  (2.27)  la  attempted  In 
the  form 

o = c®\x, t)  + e,o^^^(x,t)  -4-  £^c^^^(x,t)  ... 


u = u^®^(x,t)-*-  tu^^'(x,t)4-  g,^u^^^(x,t)  -I-  ... 


(2.29) 


After  the  details  have  been  worked  out,  the  resulting 
perturbation  functions  become 
..(0)  _ .. 


u 


(0) 


(2.30) 


i-v 

i 
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u 


(1) 


®0 

T" 

o 


.(2) 


u 


(2) 


J3) 


j^ooe(x  ooe(x  - 

u p n 

^~ooe{x  o^t)  oob(x  - o^t)j 

- rooa2x  - ooe2o^t  + 1 - ^oob2(x  o^t)-|coB2{x-OQt) J -f- 

■ o ~ -* 

Kt  [8in2(x  + Ojjt)  - Bln2(x  - c^t)j 

= ^ - ®oe2(x  - c^t)] 

- t |8ln2(x  o^t)  Bln2(x  - 

0 ^ Y+1 

where  M - - -jp  (^rT)') 

O 2 Vo 

” = ' “1“  ^2{J-i)^ 

s a^oos3(x  +c^t)  -f  a2Cos3(x  - c^t)  -f  a^00B(x  4-  o^t)  + 


,oos(x  - c t)  + a-coe(3x  -v-  c_t)  + 


agCos(3x  -c^t)  4 a,pOos(x  + ^c^t)  -i-  a0cos(x  - 3c^t)  + 
+ t|b^s1.r3(x  + c^t)  + bpSln3(x  - c^t)  4-  b^ein(x  4-  c^t)4- 
b|^Bln(x  ~ c^t)  4-  b^Bln(3x  4 c^t)  + bgsln(3x  - CQt)4. 
b^Bln(x  4 c^t)  4 bpBlnCx  - 3cQt)j  4 


4 t' 


c^cos3(x  + o^t)  4 C2Cop3(t  - c^t)  4-  3^eoa(x  4 c^^t)4 
Cji^co8(x  - Cjjt)  4 c^cos(3x  4 o^t)  4 Cgooe(3x  - c^t) 4 


o„cob(x  4 3cQt)  4 C0Bln(x  - 3c^t) 


]■ 
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where  are  constants  determinable  from  the  ocuatlone 

and  Initial  conditions,  and  where  the  function  u'-'  has  the 
same  form  with  different  constants. 

(ii) 

c'  ' has  terms  of  the  form: 

constant,  cos2x,  cos4x,  oos(2x  + 2c^t) , coe(2x  ^c^t) , 
oos(4x  ^ 2c^t) , ccs(4x  ;+ 

+ t times  sines  of  the  same  arguments, 

+ ^^2  a anna  a 

• a a B a a ^ 

The  Inadequacy  of  such  a solution  appears  as  t becomes 
Irrre,  when  the  dominant  terms  in  (2.29)  become  for  example 

c = £t(  )-^  ) -^  •••  j ■ (2.31) 

For  t = 0(g)  series  (2.31)  no  longer  converges  In  the  pertur- 
bation sense,  and  In  this  sense  the  solution  falls  since  every 
term  In  the  bracketed  expression  becomes  0(1).  The  critical 
fact  Is  that  exactly  this  order  of  t holds  at  the  region  of 
shock  development  (see  for  exairiple  page  37)  and  thus  the  usual 
perturbation  solution  cannot  hold  Into  the  region  of  shock 
development. 

2.  Perturbation  solution  In  oharaot eristic  ooordlnatee 
If  the  same  problem  (initial  conditions  (2.27))  1b 
solved  by  the  method  described  earlier  In  this  section,  the 
perturbation  functions  are  found  as, 


24^ 


u'°>  , 0 

o<°>  = e. 


= 7— v(00flfe  - 008U() 

M > 

' s ^ (008^  ■♦-OOB«i<) 

= 0) 

y for  J > 1. 

= o) 

Thus  s 2o^(A  008^-  'E  coBor) 

2 008^”  A COB«) 

are  to  be  used  In  finding 

X = x^°^-»-  Lx^^^+  ... 

t s ...  , 


u<^>  2 

JJ) 


(2.32) 


, where 


I = 

E . 


V+1 


(2.33) 


from  the  equations 

x^  s (u  +•  c)  t^ 


with 


tsOonx  = e<  = ^ 


x^  = (u  - c)  t^ 

The  rssultant  solutions  are 


S 

ot4-  ^ 
2 

- ^(0) 
°o^ 

= 

x'l> 

a 

(cOb  ^ - 008  fl(  ) A 

®o'= 

“ 

1 CM 
t 

(-008^  - c08oi)A  + (elnot-  sln^ 

x<2> 

S 

o<-  a 

(cob2«  - cob2  ^ ) (^- -y- *4" 

(2.34) 


3^(2oos«<8ln^  +2co8$Elnoi  - elnZ^  -sln2o() 
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= -^^^(coe2«x  -t-coe2p  ) (-^ + — |4^{cos  o(C08^  )2A^4- 

■— (A"  2B^  - AB)  -♦•  teri3f.  not  multiplied  by  ( - p ) , 

r ■— ^-^•(C0S3^  - 0083  «)  ( ) 3t.C.... 

r -"g— C084p  - C08<4’CX)(  ) etc... 

The  form  n8sumed  by  succeBclve  terms  is  apparent. 

In  p, articular  only  the  first  pover  of  appears,  ajid 

■ Is  the  Initial  approximation  to  t.  It  thus  seems  pos- 

sible to  continue  the  solution  beyond  the  point  whore  t becomes 
C(^).  That  this  3^  Indeed  possible  Is  shewn  In  section  2 v;here 
the  details  of  the  process  are  explored. 

At  this  point  It  is  of  some  Interest  to  remark  that 
since  (2,32)  and  (2.3^)  constitute  a correct  solution  to  the 
problem,  they  must  agree  with  the  usual  perturbation  solution 
of  (2.30)  and  also  with  the  solution  achieved  using  Riemaaia's 
function  in  all  cases  where  the  latter  two  solutions  arc  ob- 
tainable and  valid. 

For  the  case  of  the  usual  perturbation  solution  the 
agreement  may  be  demonstrated  simply  by  Introducing  the  coor- 
dinate perturbations  of  (2,33)  and  (2.3^)  Into  the  solution 
(2.38)  and  showing  that  this  causes  the  higher  order  perturba- 
tions (J  >1)  In  u and  0 to  become  zero  in  agreement  with  (2.32). 
This  has  been  checked  for  the  cases  of  u'  ''  and  . 


(2) 


.(3) 


For  the  ctse  where  the  solution  le  given  In  terms 
of  Rlemann's  function,  the  solution  may  be  shown  to  be 
(c.f,  section  82  of  [4^) 


r 


= t(  JJ(ct),  n|(^))  - 

where 


o{\)  = £g(  A )) , 


and  where  the  point 


|r  = ^(^)- 

li  r ^(o^) 


is 


the  point,  P,  in  the  (r,2)-plane  of 
the  figure  correapondlng  to  the  point  (o(  ,p)  In  the  (<v,^) -plane. 


The  functions  r and  s on  the  curve  P* 

o 


are 


r fo(^  _ 


8 (oO  = 


c_ 


8-1 


4-  6b(o()  , 


so  that  r and  e represent  the  result  of  restoring  the  ternm  of 
order  one  to  the  perturbation  functions,  r and  s,  of  equations 
(2.6b). 

In  the  expression  for  t(<y,^)  the  function  Y(A)  Is 
Rlemann's  function,  and  for  this  oroblem 


V s 


■t  s/ 


- 


where  F is  the  hyoergeoretrlc  function  and  where 

, «-l 


2?. 


Now  If,  at  this  stage,  the  prohlem  la  further  rlmpllfled 
by  asri’.nlng  that  o~\  Is  an  Integer,  the  hypergeoraetrlo  function, 
F,  beccmcr.  a finite  series.  For  exajople  by  choosing  a 1.4, 
•a^ibecomes  eoual  to  ?,  and  gives 

F - 1 - 6z  + 6z^. 

Then  using  this  F In  the  expression  for  t(oi',^)  and 
expanding  the  result  In  powers  of  E such  that 

c^t  = Ec^t^^^-t- 

one  obtains 


(0) 


o^- 


= -3/2(  ex  - p ) (cos  O'  cos  ^ ) +•  (slno^-  sln^  ) 
in  agreement  i^-lth  the  previous  solution  of  (2.34)  since  for 
8 = 1-^ 

A r = 3/2 

® = 4^^-!)  = 

Actually,  so  far  as  the  details  are  concerned,  the 
expansion  in  powers  of  ^ a.a  & perturbation  scries  which  arises 
out  of  terms  of  the  sort 

I _ I 


sCHj-tr  (o^) 


+ e(r(oO  + sl*)) 


Is  valid  only  so  long  as 

|£(r(o<)  -h  s(o^)  j - 2£K  < , 


where  K Is  an  upper  bound  on  the  functions,  r and  s 
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For  A = 


X + l 


(see  (2,8))  this  means 


e < 


Comparison  of  this  Inequality  with  the  criterion 

of  (2.18) 

’ 

shov/B  that  (2.18)  does  give  a safe  boimd  on  £ , since  for  K 1 

4>i . 

p.- 

What  this  requirement  on  £ really  means  in  terms  of 
the  allowed  Initial  depfjture  of  c from  uniformity  may  be 
found  by  solving  equations  (^.6b)  to  find 

Cq5(x)  z (r(x)  + 3(x)) 

or 

ec^lgj  ^ 05-1)  £M  » (^)Cq  . 

Thus  the  Initial  perturbation  ratio  for  ^)  = 1.4  may  be 


ec.lg 


^ )5-l  , 1 

< m “ 5 


or,  In  terms  of  the  disturbance,  to  the  density 

I 

-^  < 1.16, 

as  may  be  found  from  the  relation 
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AddltlonRl  Example 

As  further  evidence  of  the  advantegee  accruing 
to  the  method  of  coordinate  perturbation,  another  example 
will  be  considered. 

The  problem  concerns  a piston  moving  bach  and  forth 
periodically  in  a seml-lnflnlte  cylinder.  The  velocity  of  the 
pieton  Is  given  as  u^  - 6elnu>t, 

Now  If  the  usuc.l  type  of  perturbation  solution  for 
u and  c le  obtained  through  the  first  order  perturbations, 
and  If  the  flujc  of  mass  past  any  given  section  Is  Integrated 
over  an  Integral  number  of  periods  of  the  motion, 


mass  flux 


r 


e 


uat , 


(2.35) 


It  will  be  found  that  there  Is  a net  mass  flux  outwards.'  Ob- 
viously this  Is  In  error,  since  there  cannot  be  a continuous 
flow  of  matter  away  from  the  piston. 


Actually,  for  one  period,  the  mass  flux  average  Is: 
2 

= av.  density 


av.  ^u 


£ Po 

20qU> 


where 


s av.  speed  of  sound. 


(2.36) 


That  is,  the  quantity  la  of  second  order,  and  when  the  second 
order  corrections  to  u and  c(^)  have  been  obtained  and  Included, 
the  mass  flux  reduces  correctly  to  zero.  However,  this  addi- 
tional computation  required  maJces  the  task  considerably  more 
arduous.  Alec  often  In  practice  first  order  quantities  are 
used  to  estimate  say  energy  flux  and  so  on  (c.f.  a recent 


30. 


dlecueslon  and  Justification  of  this  practice  by  A.  Schook  [17]). 
and  In  this  problem  of  the  piston  the  energy  flux  cannot  be  com- 
puted correctly  from  first  order  quantities;  second  order  quan- 
tities are  required  to  make  the  energy  flux  correctly  balenoe 
the  work  done  at  the  piston. 

However,  If  the  problem  Is  described  in  terms  of  Its 
oharr.ct eristic  vrj?lables,  and  If  the  variables  x and  t are  also 
considered  In  perturbation  form  by  the  methods  of  this  Investi- 
gation, first  order  quantities  will  yield  correct  results  for 
both  mass  and  energy  flux. 

Briefly  the  results  of  the  calculations  are  as  fol- 


lows. 

Case  1)  O-dlnary  perturbation  solution 
A solution  of  the  form 

u = u^°^(x,t)  -h  iu^^- (x,t)  -«■  e u^^^(x,t)  ... 

c = c^^^(x,t)  -h  £,o^^^(x,t)  to^^^(x,t)  -f  ... 


Is  sought  for  equations  (2,28)  under  the  boundary  conditions 
that  the  velocity  of  the  fluid  at  the  piston  be  u = eslnwt, 
and  that  the  solution  contain  only  outward  travelling  waves. 
The  solution  Is: 


u<0)  . 0 

(0) 

c - 0 

u^^^  z slnu>(t  - ~) 

'^o 

(!)  ilzi6inu»(t  - 

2 Cft 


(2.37) 
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o u 
0 


°0° 


(2) 


(2) 


" -hcosa>(t  - ~)  - icos2u)(t  - ~)  +■  ^xelr.2U)( t - 


2L^ 
/ 

o 


l^) 


costUt  - ~)  ~ «cofl2w(t  - ^■:sln20){t  - •—) 


For  the  mass  flux,  from  (2.35) 
.T  T 


average  mass  flux 


= ^ r Dudc « j r £u 

“o'  V 2 


(1) 


)dt  = 


where 


= (|-)n:  ■ 

°o 


i ?P. 

2CqU> 


If  second  order  ttrirs  are  Included; 


^ J"  ( ) (u^*^U  e e,^u^^hdt  r,  0 4-  order  £- 


For  the  flux  at  energy, 
r m s Pu  s mass  flux 


given  < 


e « Internal  energy  = C^T  = |-=^  ^ (for  the  polytroplo  case) 

2.  ife 


o = 


p s pressure 


the  flux  of  energy  Is 

e)  -t-  up  = “^2^^  P o ^ « (2.38) 


The  work  done  by  the  piston  is  the  product  pu  at  the  piston. 

Then,  provided  the  second  order  terms  in  (2.37)  are 
included,  one  may  prove  that 

[average  energy  flux  at  any  point  = average  work  done  by  piston  s 

^2  ^o°o 


ige 

"1^ 


(2.39) 
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Case  11)  Coordinate  perturbation  solution 

A solution  of  the  form  (2.3)  Is  sought  for  equations 
(2.1)  under  the  conditions, 


on  the  piston  defined  as  t » of  g ; 


S 0 

f 

u<°> 

3,(k) 

- ^(1  - COBCUt) 

« 0,  k > 1 

t^^^  =0,  k>  0 

u(U 

at  the  ourve  ^ a 0 which  Is  the  characteristic  forming 

the  boundary  between  the  disturbed 
and  undisturbed  regions. 


. 0 

= ^ elncot 

= 0,  k > 1 

(2.40) 


— O V 

U S w I £- 


,(0) 


u 


u 


u 


(0) 

(1) 

(k) 


The  solution  is 
0 

slnto^  , 


.(3t) 


.(0) 


s 0,  V > 0. 


^(1)  ^ Slnco^ 


V^) 


0,  k > 1 


x^°^ 

, AO) 

<r. 

• -2"(  <x  -t  ^ 

AD 

1 . ?+i 

X 

= ui  ^ “ar 

r 

ill 

8 

(2.41) 


( 0^  - ^ ) BlnCU^  - i 0081U5H- ^ costo^ 
-(«<  - ^ ) sln(jup  - ^ COBUW  4-  ^OBUl^ 


(2.42) 

1 


The  second  order  perturbations  are  not  required. 
(Incidentally  a proof  i.'f  the  converf'ence  of  the  pertur 
batlon  series  for  this  cnee  hs.s  been  worked  out,  but  will  not 
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be  .Tlven  since  the  detells  follow  closely  those  of  the  proof 
for  the  Initial  value  problem.  Fcr  this  boundary  problem, 
given  a case  where  the  perturbation  to  the  piston  motion  Is 
bounded  by  a constant  M,  the  renulrement  on  £ le  found  to  be 

^(7  - a)*^  (2.‘0) 

For  case  11)  the  Integrations  for  mass  and  energy 
flux  must  be  carried  out  along  the  curve  In  the  -plane  which 
correctly  represents  a constant  x locatj.on.  IThen  this  Is  done, 
only  the  first  order  solutions  of  {2.^2}  need  be  used  to  predict 
zero  maos  flux  and  correct  energy  balance.  A very  definite 
saving  in  effort  Is  effected,  since  the  complexity  of  the  comjou- 
tc.tlon  greatly  Increases  at  the  higher  order  perturbations. 

Section  2:  Shock  development 

In  the  solution  to  a problem  of  one-dlraenslonal  wave 
propagation  a compression  wave  may  appear  which,  steepening  as 
it  travels,  forms  ultimately  a shock  front.  For  the  case  where 
the  shock  develops  in  a slrapls  wave  region,  Friedrichs  has  shown 
[6^  that  the  solution  nay  be  continued  with  second  order  accuracy 
as  a simple  wave  solution  through  the  shock,  so  long  as  the  shock 
Is  weak,  since  the  neglected  change  In  entropy  Is  of  third  order 
In  the  shock  strength. 

We  wish  to  show  that  a similar  approach  may  be  used 
when  a shock  forms  In  a region  which  Is  not  a simple  wave  region. 
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but  which  needs  both  sets  of  characteristic  variables  for  Its 
description.  In  this  section  the  type  of  perturbation  solution 
in  terms  of  characteristics  found  In  the  previous  section  will 
be  studied  from  the  point  of  view  of  Its  ability  to  describe  the 
flow  at  the  replon  of  shock  development  and  even  somewhat  beyond. 

In  general  the  development  of  a shock  In  a flow 
described  by  characteristic  variables  Is  characterized  by  the 
fact  that  the  mapolng  from  the  characteristic  plane  to  the 
physical  plane  ceases  to  be  single  valued;  the  Image  of  the 
characteristic  plane  folds  over  to  form  a three-sheeted  surface 
in  the  (x,t)-plane.  At  the  edges  of  the  fold  the  Jacobian, 

J = x^tp  - x^t«  (2.4t^) 

is  equal  to  zero,  and  within  the  region  the  velocity  and  density 
at  any  point  (x,t)  are  triple -valued.  This  paradox  Is  resolved 
by  the  Introduction  of  a shock,  since  the  Jump  in  the  values  of 
velocity  or  density  between  the  outermost  sheets  Is  condoned  If 
a shock  wave  Intervenes. 

General  descriptions  of  the  geometry  of  the  situation 
are  to  be  found  In  Graggs  [5],  Stocker  and  Meyer  [19] » and 
Meyer  [l5]«  The  second  problem  of  introducing  a shock  front  In  the 
correct  way  to  restore  the  right  solution  ha.s  been  studied  usually 
for  specific  cases.  As  noted  above,  Friedrichs  14]  has  considered 
the  simple  we.ve  case,  Stocker  [l8]  conelderec’  a particular 
situation  involving  a gas  of  adiabatic  exponent,  5/3*  Whltham 
[23]  » [24]  has  studied  the  development  of  shocks  in  the  spherical 
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and  axl-synimetrloel  cr.eefl  by  uilnf  coordinate  perturbations  of 
the  type  considered  in  thin  paper.  However,  he  coneldero  the 
nonlinear  deviations  of  only  one  set  of  charr.cterlstlcs , and 
does  not  consider  the  other  set,  and  thus  really  deals  with 
the  space  analop^ue  of  a simple  wave. 

If,  on  the  other  hand,  the  deviation  of  both  sets  of 
characteristics  from  straight  lines  Is  considered,  as  Is  done 
In  section  1,  considerably  moi-e  general  types  of  problcn.u  should 
be  solvab  :.e.  In  the  follov'lng  the  series  In  terras  of  charac- 
teristic variables  which  were  worked  out  for  the  case  of  a plane 
wave  In  section  1 ere  shevn  to  converge  even  In  the  region  where 
the  mapping  becomes  multiple-valued.  Then,  Just  as  for  the  case 
of  simple  waves,  the  series  may  be  used  to  extend  the  solution 
beyond  the  point  of  shock  wave  Incidence  with  an  error  denendlng 
only  on  neglect  of  entropy. 

PeveloToment  of  the  envelope  of  characteristics 

In  order  to  locate  the  region  where  the  multiple- 
covering sets  In,  the  Jacobian  (2,k4)  Is  considered.  Further- 
more only  the  cases  where  the  shock  develops  Initially  from  one 
set  of  characteristic  variables  Is  treated.  By  this  Is  meant 
simply  that  the  second  set  of  chejrr.ct eristics,  which  Is  treated 
In  the  fully  nonlinear  fashion,  does  not  form  an  envelope  at 
exactly  the  same  time  as  does  the  first  set.  Specifically  It  Is 
assumed  that  the  multiple- valuedness  occurs  along  the  ^ -charac- 
teristics whose  Images  In  the  (x, t)-plane  form  an  envelope  de- 
fined by 
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Then  from  equations  (2.1),  for  (u  > c)  ^ 0, 


(2.45) 


(2=46) 


0, 

and  the  Jacobian 

J m - X^t^  M 0,  (2.47) 

so  that  the  mapping  of  the  («,^) -plane  Into  the  (x,t) -plane  Is 
no  longer  one-to-one. 

Then,  given  a perturbation  solution  of  the  type 


c t - 0 £c  t'^^4- 


(2.46)  becomes 


■»«o' 

t ^ 


-V  £, 


2^ 


(2) 


. . 3 0 

(2.48) 


(2.48)  can  hold  only  when  the  factor,  — , of  £ becomes 

1 ^ ^ 
of  order  - so  that  the  second  term  of  the  series  can  balance 

Ci 

the  first.  To  see  what  this  means  in  terms  of  the  variables 
CK  and  ^ , consider  the  perturbation  functions  given  by  the 
general  form  (2.9).  For  the  first  terms  one  finds; 


* ! 

> 

f 

> 


X -*~R**,  [(r(0  + s(«^cl^ 


at  moat 


where  from  (2.16)  the  functions  r end  s are  bounded  by  M. 
Then 
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So  If  the  initial  functions  have  derivatives  satisfy- 
ing (2.50)  for  the  entire  range  of  their  argunients,  one  can  have 

•^c  t 

r- r — s 0 only  for 

^ 1 

(o^-g)  = 0(i),  (2.51) 


and  a shock  could  be  expected  for  lairge  oC  and/or  ^ correspond- 
ing to  large  t.  The  crucial  fact  involved  here  is  that  the 
perturbation  serlcc  converges  beyond  this  point.  If  the  power 
of  - ^ ) appearing  in  the  higher  order  terms  of  the  perturbation 
serlec  had  Increased,  the  convergence  of  the  series  under  condi- 
tion (2.51)  would  fall,  but  under  a situation  of  properly  bounded 
initial  functions,  equation  (2.18)  s'^ov:s  that  only  the  first 
power  of  (ex'  - enters  for  any  perturbation  function  and  the 
convergence  persists. 

G-eometry  of  the  shock  development  region 

The  correct  introduction  of  a shock  into  the  region 
of  multiple-valued  solutions  requires  some  exploration  of  the 
geometry  of  the  mapping  of  the  (w,^)-plane  Into  the  (x,t)-plane 
in  the  region.  In  particular  the  earliest  point  In  x and  t at 
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which  the  envelope  appears  Is  of  interest.  Such  an  extreaiuiu 
for  X is  given  by 

4*  = |Jl<K  - 0.  (2.52) 

or  along  the  envelope  defined  by  (2,^5) 


^ ^ ^ n 

dp  “ o0<  d^  “ 


(2.53) 


Since  from  the  assumption  that  a shock  does  not  form  in  the 
other  direction^ 

the  extremum  is  given  by 


^ - 
df - 


0. 


(2.56J 


Finally,  since  from  the  equations  (2.1)  for  (u  - c)  0 


the  extremum  for  both  x and  t is  given  by  (2.56). 

For  simplicity,  the  variables  of  the  problem  are  now 
shlfted  to  make  the  extremum  of  the  envelope  ocjcur  at  x s t = 0 
in  the  (x,t)-plane  and  at  0^2  ^ = 0 in  the  (o', (3) -plane.  An 
expansion  of  the  functions  x,  t,  u,  c in  terms  of  c/  and  ^ near 
this  extremum  la  then  made  in  the  form 

X = + ^11®^P'*'*20'^  ■+  ^02  p ^03J?+. .. 

t = ^10®^ 

u=  £(uoo*^"io®''*''"oiP'^*-*^ 
c = c^-l-  ...) 


(2.57) 
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where  etc,,,  are  conetanta,  and  ^-.'here  the  fact  — = a 0 

near  the  extremum,  c<  - ^ hac  teen  used, 

Alco  It  may  be  shown  th.at  at  the  extremum, 

z ~ (2.58) 

80  that  the  envelope  la  cuaped.. 

The  proof  of  (2.58)  follows  from  (2.5^>)  which  B'ows 
that  the  Iraagre  of  the  envelope  In  the  (o^,^) --alane  haa  an  ex- 
prnalon  of  the  form 

o<  = -V-0(^^)  + ...,  k = constant  (2.59) 

Then  from  (2.57) 

2 

x^  s ^ 2x^^e<  ^ P + . . » , 

or  along  (2.59) 

x^  r ^*12^ • (2.60) 

If  (2.60)  la  to  represent  x^  correctly  for  small  finite  distances 
In  ^ along  the  locus,  x^  s 0,  then  the  coefficient  of  ^ must 
vanish, 

z 0,  (2,61) 

2 3 

and  the  terms  In  ^ ^ stc.  must  balance,  e.g. 

x^2j^  ~ “33Cq^,  (2.62) 

Relations  of  the  type  (2.62)  are  demonstrated  In  a later 
example  where  the  series  (2.57)  are  given  explicitly,  (pagei4-9) . 


From  (2,61)  and  similar  conplderetlona  for  t,  the 
equation  (2.58)  is  proved  and  the  envelope  Is  cusped.  The 
series  (2.5?)  become 


X = 5^10  ^03r^  +X3QO?■^ 


t = t,.o<  ^ 


^ ~ ^10 '*'’^01  • • ] (2.57*) 

° = °0  ^ [®oo  ■*■  °10°^  ■*■  °01  ^ • 

For-  reference  In  the  dlscueslon  sketches  of  the 
mapping  near  the  cusp  are  shown  In  Fig,  1.  In  lA  Is  shown 
the  Image  of  the  envelope,  (2.59),  in  the  (^ , ^) -plane.  IB 
shows  the  envelope  In  the  (x,t) -plane  and  Includes  sketches 
of  the  triple  fold.  1C  and  ID  plot  x(^)  and  t (^)  for  con- 
stant o(  . and  show  the  Increase  In  the  area  of  the  folded 
region  as  0(  incre-'ses  from  Its  value  at  the  cusr>. 

Location  of  the  shock  wave 

The  triple  mapping  leads  one  to  examine  the  pos- 
elblD-lty  of  Introducing  r.  shock  front  vhlch  will  allow  the 
Jumps  In  vel.oclty  r.nd  density  between  the  outer  sheets  of 
the  fold  and  which  will  allow  therefore  a continuation  of 

the  solution  through  the  shocJ-c.  These  remarks  apply  of  course 
only  to  regions  of  weak  shock  since  the  entropy  change  through 
the  shock,  which  we  are  tacitly  neglecting,  affects  quantities 
of  the  third  order  In  the  shook  strength. 
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FOR  LINES  OF  CONSTANT  Q 
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The  shock  Is  to  be  located  (see  Fig,  lA  and  B)  In 
the  (x,t) -plane  on  the  upper  sheet  (region  (1))  Just  to  the 
left  of  branch  m,  corresponding  to  a location  on  the  lower 
sheet  (region  (0))  Just  to  the  right  of  branch  n.  Thus  In 
the  («,^) -plane  the  Image  of  the  shock  falls  outside  the 
region  (2)  and  behaves  near  the  origin  like 


cy  r ^ -f-  . . . , where  E > k of  equation  ( 2. 59) . ( 2. 63) 
Then  for  a point  (x,t)  on  the  shock  on  the  upper  sheet,  there 
must  correspond  the  same  point  (x,t)  on  the  lower  sheet,  and 
the  Jumps  on  velocity  and  density  between  the  two  sheets  must 
obey  the  shock  trajisitlon  equations  (Rankin e-Hugonlot  condi- 
tions [4I). 

To  this  end  let  the  shock  be  represented  parametrically 
by  a parameter  o'  , representing  distance  along  the  shock  such 
that,  In  view  of  (2.59) 

Z ') 

0/ ■ C2  <r  4-  O' ^ -4“  . . . I 

V on  the  upper  sheet,  and  (2.64a) 

^ is  d^  3“  ^2  4-  d^  <r  ^ 4*  • • • J 


o(  ss  ^ ^ • • • 

^ = ^1  ^ ^2  4 • • ■ 


^ on  the  lower  sheet. 


(2.64b) 


where  C2»  d^,  ...  are  constants. 

In  order  that  points  (x,t)  on  the  outer  sheets  may 
correspond  for  a given  , equations  (2,57')  show 


^3. 


^10' C3<r^)+  +2Tp^d^(r^+  0(<r^)  » e^tf-^)-}- 

^11  ^ ^ ■*■  ’'"  ° ' 


or,  from  the  O’  terns, 


°2  - ®2* 


(2.65) 


The  shock  trsmsltlon  eauatlonc  are  conservation 


of  mass, 


Po'^o  = = 


m. 


Tj  , U,  - U 

U = shock  velocity 


(2.66a) 


conservation  of  momentum: 

Po^o  + Pq  = ^1^1  pj  (2.66b) 

conservation  of  energy, 

2 2 

"5“  -+-  ©o  H " “5“  oT  (2.66c) 

Increase  of  entropy 

mS^  mS^  ; 3 = entropy.  (2.66d) 

In  the  solution  being  considered  entropy  change  Is  to  be  neg- 
lected which  means  that  the  last  equation  above  Is  Irrelevant, 
and  also  that  the  energy  equation  (2.66c)  reduces  to  Bernoulli's 
law  and  is  already  satisfied  by  a solution  of  the  original 
equations  ( §55,  W ) . To  satisfy  the  other  equations  use  is 
made  of  equations  (2,57*),  (2.66),  and  the  further  expansion 
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4J+. 


e = fo  + ^ 'foo  + fio~  (bi? + fo2  e *+ 


(2.67) 


and  (2.66a)  1b  evaluated  on  the  upper  and  lovrer  sheets  along 
the  shook. 


Taking  the  shock  velocity  as 
TJ  . -h  , 


( 2. 68) 


and  substituting  for  the  quantities  In  (2.66a)  evaluated  on 

the  two  sheets  of  (2.64,  a,  b)  gives  as  the  conditions  ob- 

2 

talned  by  setting  the  coefficients  of  cr  and  o*  equal  to  zero 
respectively, 


(di  - f^)  ^ fo  ^ poi'^oo)  - ° 

'^1  - ‘■i>[(tfo-£V°c”^o2  + £Voi“oi 


(2.69,a,b) 


In  these  expressions,  the  relation  d,  m f^  is  not 
allowed,  since  in  equations  (.2,64,  a,  b)  ^ must  have  differ- 
ent signs  on  the  two  outer  sheets  of  the  shock  (see  Fig.  lA) . 
Thus  one  obtains  from  (2.69a) 


XJ 


(0) 


or  since, 


''Ol  = ""Ol 


y~i  (Op 

^ °01  ’ 


u 


(0) 


(2.70) 

(2.71) 
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Then  one  may  show  from  (2.69b)  that 

^ = -^1 
TJ'^^  = 0 

iB  a satisfactory  solution  since  to  this  order  the  required 
relation  on  the  shock, 

= U Is  satisfied. 

dt/ds 


(2.72) 


(2.73) 


The  equations  for  the  shock  become 


^ = dj^  0-  + . . . 

2 I 

oL  = CgS-  -J-  . . 
^ ~ “d^(T“  -4* . . . 


on  the  upper  sheet 


on  the  lovrer  sheet, 


(2.74a, b) 


Furthermore  the  speed  of  the  shock,  that  Is  Its 
slops  In  the  (x,t)-plane  Is  the  average  of  the  slopes  of  the 
^ -characteristics  forming  the  two  bremches  of  the  envelope. 
This  Is  EhovTi  by  using  equations  (2.1)  and  (2.57')  to  obtain 


dx 


sconstant 


.-J^.u+c.c„+ 

Thus  for  equal  distances  in  ^ ( + and  - ^ ) along 


the  two  branches  of  the  envelooe. 
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r r 

= ^|^20q  + 2£,|^(Uqq+  ®oO^ 

( ^-^  ) (Uq^  + Oq^)  + o(^^)J  r 
= ®0  £(^00  ■*■  ®00^  “ 


(2.76) 


The  nomentuin  equation  (2.66b)  holde  through  order 

2 

as  one  may  show  moet  easily  by  using  the  equivalent 
relation 


Pi  - Pq  _ ^(Pj-i-  Pq) 

pi  “ fo  ~ flT  fo 


and  expressing  p and  p on  the  two  sides  of  the  shock  In 

terms  of  an  ( o(  («-),  ^((f))  expansion. 

It  Is  found  then  that  If  the  shock  wave  Is  properly 

Imposed  on  the  region  of  triple  mapping,  the  shock  transition 

eoua.tlons  (2,66)  may  be  used  to  account  for  the  Jumps  In  the 

•> 

physical  quantities.  The  equations  hold  through  order  E.<T“ 
where  Is  being  used  as  a measure  of  shock  strength  since 

for  example  from  equation  (2.6?) 

po  foip"'’ 

which  along  the  shock  becomes 

f = fo  + s[Poo  °<»^>] 

shows  that  a Jump  in  density  is  of  the  order  of  a Jump  In  £.0*. 

As  far  as  the  error  due  to  neglect  of  entropy  change 
Is  concerned,  one  would  estimate  It  as  the  cube  of  the  shock 
strength,  or  (£<T)^,  but  actually  the  error  Involved  Is  larger 


since  the  ehoch  transition  equetlone  sre  satisfied  only 
through  order  £<r  and  Involve  an  error  of  order  S.<r , 

To  find  the  actual  error  Incurred  by  pushing  the  solution 
a distance  x beyond  the  point  of  shock  formation,  equations 
(2,57')  and  (2.74a,b)  are  used  to  show  that  since 


X — < ^2  g~Z  ^ 

a 

a:'  error  of  order  corresponds  to  an 


error  6x  (2 

111  traveling  a distance  x along  the  shock. 

Illustrative  example 

The  example  used  In  the  first  section  1l!  reexamined 
here  from  the  standpoint  of  shock  development.  The  periodic 
initial  distribution  of  the  problem  shows  th.at  not  one,  but  an 
infinite  set  of  shocks  will  develop  in  a repeating  pattern 
from  the  Infinite  Initial  set  of  compression  waves.  Further- 
more envelopes  will  occur  for  both  the  of  - and  p -sets  of 
characteristics  and  a sketch  of  the  pattern  of  envelopes  is 
given  on  page  51.  However,  since  this  study  does  not  purport 
to  study  Interactions,  a single  developing  shock  wave  will  be 
Isolated  and  used  to  Illustrate  the  points  made  earlier  In 
this  section. 

From  equations  (2.34),  the  envelope,  = 0,  la 
found  to  be 


A 


£Asln^  aln^ 


(2 
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vhere 

m (2i+’4B)  coepoo8«<»  - 3fislnA8ln^  - 4-  ®)  cos2o< -t- 

(-  I “ ^2s)ooa2^  ■>■  (a^.|)4-  i2A±Il,.  fcjlrj^d-*- coe2p). 

2A 

If  the  problem  ie  made  more  definite  by  ohooeirig 
y = 1.4 
A = 3/2 
g - 1 

r _ T /T  nn 

^ S w 

then  the  envelope  becomee 


„ ^ 6+^20_  + coea  -I-  .J.  Q(g2, 

' 3Bln^  sin^  8ln^ 


(2.79) 


where 


a 7coepcoso/  - 3slno^8ln^  - cob2<*  - 4cos2|9  + 24- -^8ln|5(l-t-coB2^ 


An 

The  cusps  of  the  envelope  are  located  at  0 and  choosing  the 
the  branch  which  has  a cusp  near  ^ a ‘^/2,  the  cusp  is  at 

= b9»2 
:tc-  3‘>.9 


Then  if  the  problem 
the  origin  In  both  the  (of,p)- 
expi.nBlons  of  the  type  (2.59) 
The  equation  of  the 


33.5 

is  shifted  to  put  the  cusp  at 
and  (x, t) -planes,  the  following 
and  (2.57')  apply, 
envelope  In  the  (Q^,^) -plane  near 


the  cusp  has  the  expanelon. 


Lio 

y * 


3i.9(^  ...  . (2.eo) 

Tre  quantities  x,  c^^t,  u,  and  c near  the  cusp  exhibit  the 
following:  behavior. 

X c^£  ... 

c^t  cs  +.52<«  + .008«p- .OSp^-h  ..  . (2.81) 

+ .06o^  + .005^ -h./f99of  ^'-h  ...1 

f-C:'  1 -h  1^0.988  - ^ - .06«  + . 005^  ^ - .499«  ^4-  . . !j  . 


Th^t  Is,  the  expressions  (2.81)  do  astnine  the  form  predicted  In 
(2.57') » and  In  fact  relations  of  the  type  (2.62)  are  verified 
since  In 

Xlik  = -3xq^ 

The  values  are 

k r 31.9 

= -.008 

Xq^  s +.08 

and  the  equality  holds  to  within  the  accuracy  of  the  calcula- 
tion. 


The  error  Involved  In  this  particular  problem  may 
be  estimated  again  from  the  failure  of  the  shock  transition 

p 

equations  to  hold  through  order  £.  <T " . From  the  first  of 
equations  (2.81) 


£!T 


J5S  16^^  +-0(  ^^)  as  16(t2  + O(cT^) 

3/2  £ 

JS  -'  ^ W 


so  that 


(2.82) 


and  the  solution  may  be  continued  for  a dletanoe  x along  the 
shock  with  an  error  given  by  (2,82), 

As  etated  previously,  the  actual  flow  pattern  In 
the  (x,t) -plane  Is  a configuration  of  Intersecting  compression 
and  expansion  waves,  so  that  a.ependlng  on  the  tim©  of  shock 
development  the  compression  wave  traveling  to  the  right  and 
steepening  has  encountered  a certain  number  of  leftward- 
traveling compression  and  rarefaction  waves.  Prom  the  fact 
that  a shock  develops s at  s time,  t,  of  order  ^ one  sees 
that  the  patterns  of  interaction  vary  with  £,  and  that  in 
fact  for  certain  epsilons  the  shocks  traveling  in  different 
directions  would  tend  to  form  at  the  same  time,  that  is,  the 
envelopes , 

-+  0(  £ (2.83) 

sln^ 

coap  ^ ^^Eeoaor  , 

~=0  : ~ + A h --^4-  o(  £ ^)  (2.8^0 

c Asln«?  sinc*  slnof 

wliere  R^(«/,^)  r ap{^,«=<) 

would  each  have  a cusp  at  the  same  values  of  and  ^ . In 
Fig,  2 on  page  ^1  are  shown  sketches  of  the  pattern  In  the 

(x, t)-plane  and  of  the  effect  of  £ on  the  pattern  of  Inter- 
action. The  Initial  distribution  of  c Is  superimposed  on 
the  dlagraiCB.  For  the  problem  discussed  itove  with  £ - 1/100, 
the  shock  development  is  delayed  so  long  that  many  Interactions 


^x 


0< 


OOSC<  '»• 


£5ein^ 


Binp 
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have  oocurred;  for  6'^*^  on  the  other  hand  only  one  Interaction 
has  occurred,  and  for  some  t ^ 1/6  the  pattern  of  Fig.  2c  oocara. 
This  last  estimate  cannot  be  made  with  assurance  from  the 
equations,  since  the  £ Is  larger  than  that  allowed  on  page  28 
from  convergence  considerations,  but  obviously  there  will  be 
a first  £,  to  give  such  a region,  and  there  will  also  be  a 
sot  of  Increasing  £*s  Involving  such  regions.  The  behavior 
of  the  solution  In  these  neighborhoods  has  not  been  explored 
In  detail  and  might  prove  of  considerable  Interest. 

B.  Cylindrical  and  Spherical  Wave  Propagation 
Introduction 

The  success  with  which  the  coordlnate-pertiirbatlon 
type  of  solution  was  able  to  describe  the  region  of  shock  forma- 
tion for  the  case  of  the  plane  wave,  leads  one  to  hope  that 
this  type  of  solution  might  bs  extended  to  describe  shock 
regions  for  the  spherical  and  axi-eymmotrlcal  cases.  However, 
in  the  olane  wave  case  the  right  to  carry  the  cclutlon  ur>  to 
and  through  the  shock  depended  on  a definite  proof  of  the  con- 
vergence of  the  series,  and  unfortunately  such  a convergence 
proof  has  not  yet  been  found  for  the  higher  dimension  cases. 

Thus  the  right  to  extend  the  solution  for  these  cases  must 
be  held  In  abeyance. 

Nevertheless  It  seems  probable  If  not  yet  Justified 
that  the  extension  can  be  made,  and  In  fact  G.  B.  Whltham  (23)  %[24l 


c 
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by  allov/inp  a t>erturbatlon  to  the  direction  of  one  eet  of 
chr.r?ote..'lstlc6  haa  had  considerable  success  In  studying  the 
regions  of  shocks. 

Pone  aspects  of  the  matter  wli.l  be  considered  here. 

In  section  1 the  general  form  assiuned  by  tlie  equf.tlons  for  the 
sequence  of  perturbation  functions  will  be  given  together  with 
a brief  discussion  of  the  type  of  convergence  conjectured  for 
the  series  solutions. 

Section  2 compares  the  results  of  section  1 with 
those  found  by  n'hltham  In  [23]  . 

Section  1;  Perturbation  solutions 

In  terms  of  characteristic  variables  e<  and  ^ , the 
equations  describing  spherical  or  cyllrdrlcal  wave  motion  be- 
come ( % 23  ,(^] ) 

= (u  ■+  c)t^  (-J^)  [ u^  -f  t.,1  + c^  - 0 

(2.85) 

= (u  - c)tp  [u^  - t^j-  Cp  = 0, 


where  n Is  the  number  of  dimensions  Involved,  and  where  x and  t 
are  the  space  and  time  coordinates  respectively,  and  u and  c 
the  velocity  and  speed  of  sound. 

For  simplicity  the  initial  conditions  will  be  given 
for  a gas  Initially  at  rest,  fuidas  before  the  Initial  line 
will  be  defined  by  setting  both  cv  and  ^ equal  to  x. 

Then  cno(=p  , x s o(  sa  ^ 


t B 0 

u ■ C 


(2.86) 


+ £.  c^^^  (x) . 


c = 


If  a perturbation  solution  Is  sought  In  the  form 


t;  — €.2^{2)-h  ••• 

u . u^°^+  £u*'^^4-  ... 

e « 0^°^  + 4 + . . . 

frorc  (2.85)  and  (2,86)  the  terns  free  from  L are  found  to  be 
^(0)  ^ 


(2.87) 


0 

o"  “2 
u^®^  « 0 
- Cq. 


(2.88) 


Using  (2,88)  the  equations  for  the  first  order  per- 
turbation furxctlons  u^^^  ijid  become 


(2.89) 


(2.90) 


At  this  stage,  both  to  expedite  the  calculations,  and 
to  obtialn  a solution  In  a form  comparable  with  vnilthan's  work 
on  spherical  blast,  the  problem  v'lll  be  cpanlellzed  to  the 
spherical  case. 


Then  (2.90)  becomes 


(1)  ^ 


cX-\~  8 


(o<+f») 


= 0 


which  nay  be  shown  to  liave  the  eolutlon 


(1)  _ 2[f  (of)  H-  g(p)] 

^ “ of+-  ^ (e»f-*-p)2  ’ 

v:here  f and  g are  arbitrary  functions  determined  from  the 
Initial  conditions.  Finally  using  the  Initial  conditions 
and  the  equations,  a conjt>lete  set  of  solutions  for  the  first 
order  uerturbations  Is, 


u 


(1) 


1 


df  (3<) 

. 1 

2 f. 

dof 

dp  J 

(o<+p)2  L' 

Lil  _1_ 

\^(P)  ^ 

df  (o/)  1 

2 

“2^J 

[f(c<)  - f(p)J 


(1) 


>0^ 

4 


(g-J) 


OV>^ 


M - Ife)  _ iife} i, 

if  Zot  4b  ^ 


— ^y) * ®o  that  the  equation  for 

hecones  c^^^  - (o<) -h  (p)j /(<xtp) . 

Hov/ever,  at  the  next  step,  when  the  higher  order 
(2) 

perturbation  functions,  u ' etc.,  are  to  be  found,  the  tern 


where  — 


2.91) 


2.92) 


{ 


(2.93) 


(n-l)uoti 


In  (2.85)  IntroducwB  non-homcgeneous  terms  Into 


the  equations  and  mates  the  solutions  of  the  equation  much 
more  difficult  to  obtain.  It  is  exactly  this  difficulty  in 
finding  the  expressions  for  the  solutions  which  has  deterred 
the  proof  of  the  convergence  of  the  perturbation  series. 

At  any  rate  the  form  of  the  equations  for  the  higher 
order  perturbations  can  be  given  for  a general  n“**  order  per- 
turbation. In  view  of  the  fact  that 

1 1 li  . 

I = 7^1^  • ^ '•  V ‘ J ^ 


i , ••  • >0 

and  with  x^^^  s , the  general  equations  (2.85)  become, 


[“i"' 


)Czl 


u. 


,(n)  ]^in) 
j 


V. 


zit -1)  \ 


*+  p ^ 

o+q+r-^J  _ n 

P,q.r,j  < n 


0» 


hiv«es 


(n) 


\ * I r — 

-2-u<">!  - c<“>  - 


p+q4-r+J  3 n 
p,q,r,J  < n 
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Now  even  though  a proof  of  the  convergence  of  the 
perturbation  series  has  not  been  found,  It  seems  worthwhile  to 
diocuss  the  region  where  a shock  might  be  expected  to  develop, 
and  to  attempt  an  estimate  of  the  size  of  the  terms  In  the 
perturbation  series  in  such  a vicinity.  To  this  end,  as  In 
the  plane  wave  esse.  It  wl?^l  be  resumed  that  only  one  set  of 
characteristics  folds  Into  a triple  mao-ilng  at  one  time.  The 
region  will  be  considered  asymptotically  In  terms  of  only  one 
characteristic  varla.ble. 

Pli  st  of  all  if  the  perturbation,  to  the 

quantity,  c.  Is  to  vanish  initially  for  large  distances,  (2.93) 
shows  that  asymptotically 


— -w-  constsnt, 

dy 

.(1) 


(2.95) 


Then  from  the  eouatlon  for  x'  ^ In  (2.92),  under  the  assumption 

(2.96) 


^ >> 

one  finds  . 

Considering  therefore  the  asymptotic  form  for  the  expression 


X = £x^^^- 


0(p)  -h  0(  6.  io  ^ ) -t* 


• • • > 


(2.97) 
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W©  note  that  for  the  deviation  from  the  linear  case  to  be 
great  enough  to  allow  the  charge terletlce  to  croee,  the 
second  term  must  be  of  order  1,  or 

(2.98) 

(c.f,  a similar  conclusion  by  Whitham,  h^{23\  ).  The  question 
then  arises  as  to  the  order  of  the  higher  perturbation  terms 
In  this  vicinity;  that  Is,  does  the  series  converge  In  such 
a region?  In  answer  only  the  unproved  conjecture  caJi  be 
stated  ci^ y uip oo b 1 caxxy 


x<")  c t'"'  ~ 
o 


n > 1. 


(2.99) 


The  heuristic  argiment  leading  to  the  relations  (2.99)  is 
based  on  a consideration  of  the  non-homogeneous  terms  In  (2.9^). 
For  example  the  equation  for  the  second  order  perturbation  has 
terms  of  the  form, 

= Jl)„(l)  ^ 


u^-^^  + 


(<y+p) 


u 


and  If  asymptotically 

^(l) 

\ • 

then  asymptotically  +■  , , . ^ , 

and  then  from  this  without  being  able  to  Justify  the  step,  we 
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eesume 


(2.100) 


Relylnr  on  tho  conjectured  relrtione  (2.99),  the 
perturbation  series  would  beve  the  asymptotic  orders. 


ft  2 

X 0(p)  4-  0(£^f(^)  4-  0(t^  ^^)  -f  0(g.3-^J^)  + ...,  (2.101) 

or  In  the  vicinity  Indicative  of  possible  shock  devel- 

opment , 


X O(^)  + 0(1) 


0(p 


0(-> 


) +- 


(2.102) 


giving  a series  stil],  convergent  In  a perturbation  sense. 

The  work  done  by  VTiltham  nontp.lns  perturbation  series 
whose  behavior  at  large  distances  agrees  with  tho.t  found  in  this 
section,  end  Is  discuseed  In  the  following  section. 

Section  2;  Comparison  with  Whltham's  work 

Whltham,  In  his  paper  on  "The  Propagation  of 
Spherical  Blast,"  [S3] , makes  the  point  that  in  a strictly 
linearized  wave  propagation  problem  the  characteristics  are 
straight  and  parallel,  and  hence  cannot  cross  as  they  would 
in  actuality  do  near  regions  of  shock  development.  To  correct 
the  situation  he  proposes  that  the  solutions  for  velocity  and 
density  be  considered  as  functions  not  of  the  linearised 
charr.ct eristic,  c^t  - x,  but  of  a true  charaoterlstlc  varl- 
able,  z,  where  z Is  determined  to  be  truly  constant  along  a 
characteristic.  The  IjnRnr  characteristic  la  to  be  cor- 
rected by  the  expansion. 
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o^t  - X - zlogx  - h(?  ) - (m^(z)logx  + +■...» 

equation  (9)»  1^3]  , and  this  expansion  together  with  expan- 
sions of  the  type, 

velocity,  u s -c^f(2)x^^  + j^(b^(z)logx  +b2(z))x"^  + 

(c^(z)logx  4- 02(z))x“^  + ... 

are  to  be  used  to  find  the  solution  of  the  probleir. 

In  other  words,  a perturbation  to  the  direction  of 
only  one  set  of  chf racterlstlcs  is  used  and  it  is  aBSuriied 
that  there  is  no  significant  nonlinsarlty  In  the  other  set. 

Such  an  approach  differs  of  course  from  that  v;e  have  used 
earlier  In  this  report  where  perturbation  series  were  con- 
sidered for  both  variables  x and  t In  terms  of  both  sets  of 
characteristics.  In  our  ease  a correction  is  automatically 
effected  for  both  characteristic  directions,  but  for  the  spherical 
case  the  problem  tends  to  become  more  complex  than  for  the  plane 
wave  propagation. 

The  solutions  which  Whltham  finds  using  the  expansions 

above  are 


. _(C^)  , Co 


-'o 

c(^ualton  C'i) 

_ / ^ is 4-...  , 

X 


X 
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where  and  B2  are  arbitrary  oonstante,  k = 

kg  = ^ and  where  the  Tunction  h(x)  Is  determined  from 

Initial  oonditions. 

On  the  other  rand,  the  Bolutions  for  theae  functions 
whicii  we  find  using  our  kind  of  perturbation  approach  are 
equations  of  the  type  (2.92),  which,  before  enforcing  Initial 
conditions,  become 


u 


(1) 


d<x  Ie> 


0/  ^ 


(o  + p)2 


(2.103) 


10g(«*4(6)  [g((S)+  f(^)j 


(i^y) 


J 


(2.104) 


where  in  a function  to  be  determined  from  the  conditions 

of  the  oroblem. 

Now  in  order  to  be  able  to  compare  ^fhlthani’s  results 
with  ours,  several  clarifications  must  be  made.  Primarily  of 
course  since  his  results  are  based  on  outgoing  waves,  our 
result  must  be  simplified  to  contain  functions  of  the  charac- 
teristic & only.  Further  since  we  have  results  only  through 
the  first  order,  the  comparison  can  be  made  only  at  this  level. 
The  equations  to  be  compared  become 


u s 


X - c^t  = 


dg(^) 

a 6 _ 2g(^) 


(o<4-^) 


2 


+ 0(  £^) 


(2.105) 


£ 


t^) 

(2.106) 
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IfhithAn;  dose  not  expressly  state  that  his  solution 
Is  a perturbation  one,  but  the  fact  that  he  oenslders  a derla- 
tlon  of  the  characteristics  froa  those  glren  by  the  linear 
solution  Implies  a perturbation  technique.  Thus  from  his 
equation  (13) » the  quantity  z vhloh  is  constant  along  a character- 
istic also  must  act  as  a small  order  perturbation  quantity 
(see  also  the  discussion  In  section  4 of  his  paper)  and  so  we 
expand  z In  a perturbation  series 


_ c 7.  ^ ^)  f ■ 

- - 'T' 


.2. 


r(2)  ^ 1 

< \ p/T  * 


v<c.  xu/j 


From  the  same  sort  of  oonslderetlons  his  terms  of  order  one 
In  both  (c^t  - x)  and  u are  set  equal  to  zero.  (Reetrictlon 
of  the  problem  to  one  of  zero  Initial  velocity  does  not  alter 
the  comparison.)  Prom  the  terms  of  order  one  In  his  equations 
then 


B,logx  4-  , 0 

-b 

B 

Bj^logx  - - k j T?h'(!;)di$  - Bg  s 0, 

0 


(2.108) 


BO  tloat  B^  s 0 and  B2  must  be  such  that 

■^h«(^)d^  . -B2-^  (^)  + ...  (2.109) 

u 

In  fact  it  Is  true  that  becomes  zero  In  the  problems 
V71tham  considers  and  that  (2.109)  or  an  analogue  thereof 
appears  in  his  section  3 In  the  form 
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= Js: 


!'>,2 
V - 


zh(  z) 


s -l/2kb' 


Using  equations  (2.107),  (2.109),  end  also  the 


expansions 


X = 4-  ,p)  4-  O(f^) 

h(z)  = H^°^(f)+  £H^^^(^)^-  fc^’(^) 

In  his  equations  (12)  and  (13),  gives: 

2c  kz|j!  4c  v1!^^^0)  , 3''  2 

— 2 — j£jL 2 L_  ^ j.n(  c^\ 


«t<4-  p 


♦p) 


(2.110) 


( o T n N 


t-  C,i 


(2.112) 


Then  upon  comparing  these  results  with  ours.  It  Is 
found  that  (2.111)  agrees  with  (2.IO5)  and  (2.112)  agrees  with 
(2.105)  provided 


H^(z)  = ^ 

= SEj^Cp)  + log  2 

20^^^^^)  r g(^) 


This  is  a consistent  set  of 
ing  (2.109)  one  finds 

lczh‘(z) 


relatione  since  by  dlfferentlat- 


-¥ 


(2.1lJi) 


64. 


and  then  using  (2.110)  and  (2.113) 


dt  _ 1 

^ h*(z) 


t 


80  that 
kzh-*(  z) 


)h’(z)  = k£ 


and 


and  (2.114)  holds  under  the  conditions  of  (2.113). 

Thus  the  first  order  solutions  are  In  a^eeu'^it, 
and  .^Itbourh  ve  do  not  have  higher  order  solutions  available 
for  ccnparlson,  it  Is  possible  to  show  that  Whlthrin's  eolations 
asymptotically  take  on  the  form  (2.101)  estimated  In  section  1. 
Asymptotically  his  equations  (12)  and  (13)  become 

u zO(i)  -f-  ^g-^)  + . . . 

c.  t ~ zOdog  x)  + z^0(^°§  ")  -f-  . . . 


in  good  agreement  vlth  our  estimates 

u ->tO(i)  £2o(12£1^) 

V f ^ 


H“  . • . 


O(^)  *>-  fOdog^  ) -V  ^ 


Evidently  a good  deal  of  work  remains  to  be  done  In  the  fopll- 
cation  of  coordinate  perturbation  to  spherical  (and  cylindrical) 
wave  propp^ation  problems.  Not  only  would  general  convergence 


proofs  be  desirable  for  a process  Involving  both  sets  of 
cliers.cterlstlcs , but  even  an  Investigation  of  the  form 
of  the  next  order  pcrturbptlon  solutions  vrould  be  usef  il. 
Partlculer  yi’cT.' lenc,  ccuic  be  studied  v/ltt  special  referene 
to  those  Involving  the  develormert  of  shock  waves  so  thr.t 
the  perturbation  solution  could  be  compared  with  the  solu- 
tion obtained  by  other  means  (e.g.,  noiterlcel)  both  In  the 
ordinary  regions  of  the  ilc\-  and  In  the  neighborhood  of 
the  shock. 
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THE  ELUPTIG  CASE 


Introduction 

The  applloftblUty  of  the  coordinate  perturbation 
t3roe  of  solution  to  problems  of  an  elliptic  nature  has  been 
studied  by  Investigating  the  problem  of  the  flow  pest  a thin 
airfoil.  In  the  usual  thin  airfoil  theory  the  solution  ob- 
tained tends  to  bree.:  down  near  the  leading  edge  of  the  air- 
foil and  30  In  particular  one  would  wish  to  correct  the  solu- 
tion In  this  neighborhood.  Llghthlll  considered  the  problem 
for  the  Incompressible  case  [111  and  found  that  the  first 
order  solution  could  be  rendered  uniformly  valid  by  Imparting 
a constant  shift  to  the  coordinate  system.  Furthermore  he 
iciplled  that  for  the  higher  order  solutions  any  tendency  for 
slng’alarltles  to  arise  could  be  suppressed  by  finding  suit- 
able higher  order  perturbation  functions  for  the  coordinate 
system. 

In  th-ls  chapter,  however,  It  Is  shown  that  the  per- 
turbation I'orctlons  required  to  render  the  solution  valid  at 
each  stage  have  Increasing  orders  of  singularity  near  the  lead- 
ing edge  of  the  airfoil  so  that  they  do  not  represent  a satis- 
factory coordinate  transformation.  Thus  although  p.  one-stage 
correction  does  effect  an  Improvement  in  Lhe  unllorm  validity 
of  the  velocity  fields  for  the  Incompressible  case,  the  general 
appjdcablllty  of  the  method  st  the  next  stages  cannot  be  af- 


firmed. 
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As  far  p.t  the  compresstble  case  Is  concerned,  even 
the  one-cta^e  correction  proves  unsatisfactory,  and  so  In 
Ceneral  it  dees  not  seem  possible  at  this  point  to  recommend 
the  method  of  coorciinate  poi-turbatlon  as  n pantcea  for  prob- 
lems of  an  elliptic  ciiaracter. 

Section  1;  Inoompresslbe  case;  Difficulties  with  Ll.^hthlll*  s 
solution. 

In  the  following  it  will  be  shown  that  if  a coor- 
dinate perturbation  of  the  t^pje 

X = X-b  (X,Y)-*-  (X,Y)-v  (X,Y)-b  . 

. _ V (3.1) 


la  to  be  used  to  correct  the  perturbation  solution  to  the 
problem  of  Ineoinpressible  flow  past  a body  whose  leading 
edge  at  X = 0 behaves  like 


y 


(3.2) 


Then  the  process  will  be  unsuccessful  foi'  3 1.  It  will 

happen  that  near  x • 0 the  pertiirbatlcn  functions,  have 


the  behavior; 


2s-l 
•«-  » 


9 


(3.3) 


£,o  that  all  the  functlonfl  will  be  positive  powers  of  x as 


X — ► 0,  only  for 

,?im  8 > or  s » 1. 

n-^«* 

Llghthlll  [llj  took  the  case  s r V2  so  that  his  x^^^  1b  a 
constant,  but  his  higher  order  perturbation  functlona  become 
singular  and  the  process  breaks  down. 
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In  a general  application  of  a coordinate  transformation 

of  tbe  type  (3.1).  the  perturbation  functions  for  the  y-coordlnate 

would  have  to  be  restored,  and  the  differential  equations  and 

boundary  conditions  on  both  sets  of  functions  and  y^^^ 

determined.  In  this  particular  problem  It  may  be  shown  that 

the  equations  can  be  satisfied  by  setting  the  pair  and 

y^^^  to  be  harmonic  conjugate  functions,  but  that  It  Is  also 

possible  to  set  y s Y as  Llghthlll  does,  provided  that  higher 

order  perturbations  In  the  velocities  (potential)  are  retained. 

The  matter  of  greater  Interest  la  the  boundary  condition,  and 

(2) 

In  fact  Llghthlll  determines  the  function  x solely  on  this 
basis. 

In  his  equation  (28),  [l3^  he  sets  u(x, y) u(X, y) 
and  v(x,y)  — v(X,y),  where  u ard  v are  the  velocity  components 
In  the  X and  y directions  respective?.;/,  smd  then  takes 


hy 


F'(X  e^x^^^)  on  the  airfoil  defined 
£F(X+ 


where  for  simplicity  we  hpvc  taken  the  camber  equal  to  zero, 
and  have  further  neglected  the  function  agreeing  with 

Ligl'itlilll  that  It  is  not  necessary.  The  point  of  the  step  Is 
to  Improve  the  usual  perturbation  solution  (u(x,y),  v(x.y)). 
which  Is  not  valid  uniformly,  by  e.xpressing  the  solution  in 
terns  of  new  coordinates  In  an  X-y  plane  and  then  by  deter- 
mining the  appropriate  coordinate  transformation  between  this 
plane  and  the  physical  plane  which  will  send  the  approximate 


solution  Into  a uniformly  valid  one.  IJov  in  view  of  this 
fflapplng'  idea  It  oeems  tippler  to  express  the  problem  In  terns 
of  a stre.aiE  function  since  in  the  first  place  the  liaag-e  solu- 
tion in  the  (x,y) -plane  can  be  clearly  visualized,  .and  secondly, 
only  oiifc  dependent  perturbation  solution  is  required.  Using 
the  otreaS.  funetiou  It  is  possible  to  'i-iupiicate  the  steps 
and  30tiv,atlon  of  Llghthlll's  pr.per  and  to  shov  how  difficul- 
ties. occur  at  the  next  order  perturbation  which  he  did  not 
consider. 


For  the  stream  function, 
the  boundary  conditions  become 


where  u 


y 

u 


I y r fcF(x) 
I 0 ^ X ^ c 


-V 


X* 


CxssO)  a leading  edge 
Cxsrc)  - trailing  e^dge 

or  ''i'lXjy'  a const.arit  s 0. 

The  aitalo.gue  of  Llghthlll’s  equation  (28)  is  then 

^(X,y)  = constant  = 0 on  y = £F(x)  = eF(X -»• 

Kow  for 

y +t  y^^(X.y)-»-  (X.y)-b  ... 


(3.^1 

(3.5) 


and  for 


"t^^^(X.y)  s + y 


and  for 


^ys)(x,0)  ^2  ^2^(1)  ^ 

^ 2 -j  2 

^ y s:niall. 


by 


y -r  = €F(X)-»-  £^x^^^F*(X) 


one  finds,  as  he  did,  a set  of  bo^uralary  equa.tior!«  of  the 
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following  t^ee,  derived  from  setting  the  coefficients  of 
powers  of  £ equal  to  zero. 

Boundary  conditions 


Tsrss  order  E 

2 

Terms  order  £ 

3 

Terms  order  t 

4 

Terms  ox^ier  % 


(X,0) 


-F(x) 


0 ^ X ^ c 


-F(x) 


-S' 


-F- 


■»y 

■©y 


"ay 


i X,0 

3^.(2)p, 

3^(3)p, 


2 -Sy2 


p2  -32^(2)  p3-j3l(/(l) 

2 -»y^  “3)  -»y3 


In  these  equations  the  arbitrary  functions,  yS^\ 

etc.  are  to  be  chosen  to  suppress  Increasing  orders  of  slngu- 
iejrlty  (nesr  x ■ 0)  appearing  In  the  successive  'B.  It  is 

necessary  to  estimate  the  order  of  the  terms  Involved  and  we 
do  this  for  a general  case 

^{Y\  r vB/ri  ,ji.  w V ^ ^ a r\ 

r va;  — jx  v»Q  ' * ° ' -• 

Prom  the  terms  of  order 

V^^(3C,0)  = -P(X)  a 'X®(Fo  + FjlX  4.  ..  .), 
so  that  the  completed  complex  potential  function, 

3(1)  = vp(i) 

which  Is  to  be  used  for  finding  — In  this  neighborhood 

V J' 

behaves  like 

= -1FqZ®4“  ...  , Z = X + ly. 
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Then  frora 


~SZ~  = -tj  ^ ^ 


-»  X 


-iFqb: 


.B-l 


Is  regiilE-x  In  the  nclg'hborhaoi  (as  y— ► O, 


one  finds  _ 

X-**  0) . Slulleo^lY 

,2p(l)  -i2Tt/(i)  -i2Ti,(l)  _ - „ , 

i^  - - ■ » — . -IFqsCs-DZ®'*-  + 0(2®”^) 


dZ* 


■^X^y 


or 


ilSD 
y— ^ 0 


^24.(1)  ^^2-(l) 


■ay" 


“a  x" 


P^eCs-DX®”^-!-  . 


and  in  ttim 

. ->3^(1)  « o 

ilM  — = — r 0(X®"^), 

y-^^O  Dy-^ 


Using:  these  eTalnationa  In  eauLations  (3-6)  shows 
that  whereas  le  finite  as  well  on  the  boundary 


<a 

(both  X ) , an  Increase  In  order  of  slng»ilarlty  appears  la 

3 


the  £ equation  date  to  the  term: 

'o 


p2^2^(l) 


^ — • FqS(s-1)X 


(2) 

A is  to  bs  chosen  to  cancel  out  this  hi^^hex*  order  singularity, 


or 


2 ■»_2 


.(2) 


Fo^(e-l)  ^8-1 


2 


£ *- 


Then  froa  the  terms  of  the  £ equation  (3*6), 

-s»3y(13  -^2^iz)  ^2^{l) 

since  — Is  not  sore  elagular  tl-jan  -g — ^ 5 — , 

y-^  "3  y"  "3  y- 

tlie  meet  singular  terns  Involve 

_(3)^=  y2-^2y(2) 

“Jt  jr  * - * = — - 


2 -it 


T-  r 0 


giving 


:<3)  _ q(3^8-1j^  which  Is  not 


more 


.(2) 


singular  than  z 

However,  the  equation  of  tersms  of  order  g ^ will 
contain  a tern: 


«y 


to 


(Ijf) 

be  balanced  by  a term,  -x  F*  so  that 


xC»)  . oix^’-h 


In  aoHsary  then  for 

X = X-**  f-M^^(X,y)+  A^^^(X.y)  4-  ... 


as  y- 
X 


0 

O'* 


A2) 


j2b-1 


x'3>  . 0(x2“-l) 

X'***  I od"*'"’) 

,(5)  X od'*'-’) 


x(6)  ^ 0(X^“-5) 


for  S - 1/2 


-2 

.(2)  lO 


x^^^  s 0(1) 

= 0(|) 

= 0(J> 

x<°^  = O(^) 


(3.?) 
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^Thun  TAg-hthlll's  determination  that  for  8 s 1/?, 

(?)  ^0^ 

x'  " ■ quite  coi-rect  and  the  correction  servee  to 

r-educf  the  error  In  velocity  at  the  lealir-T;  edge  from  0(1) 
to  0(  f X*  Hov^ever  unfortunately  the  usupl  perturbation  scheme 
Is  such  that  each  higher  order  per-urbation  function.  "V 
Introduces  a higher  order  singularity  In  thle  region  (c.f. 
equatlcnc  (3.6)),  and  the  coordinate  perturbation  functions, 
cannot  suppress  the  singularity  v/ithout  carrying  It 
themselves.  The  result  Is  th-.t  near  the  leading  edge  the 
perturbation  series, 

X = X + £^x^^^(X,0)+  t^x^^^(X,0)  -»■ 

p 

evaluated  say  for  X « 0(  ) , (which  holds  at  the  Ir-age  in 

the  (X,y)-plane  of  the  leading  edge),  becomes  (for  s = 1/2), 

X rr  0(  £ ^)-^  0(  £^)  -^0(  £^)-^  0(  £“)-*'  0(  £ ^)  + 0(  £“)  -+  ... 


ar.5  the  r^erturbatlor.  Is  no  lonTcr  valid. 

Alno  for  s ^ 1/2  from  (3. 7),  the  perturbation  func- 
tions will  be  positive  powers  of  X only  for  s ^ 1,  that  Is 
for  profiles  with  a finite  (or  zero)  slope  at  the  origin, 
and  so  a convergent  scheme  can  be  evolved  only  for  special 
simple  cases  where  the  leading  edge  Is  not  blunt.  Thus  at  this 
stage  of  the  Investigation  no  way  seems  avt.ilable  for  extend- 
ing the  method  to  a i?:enerally  useful  scheme. 


♦ The  Improvement  In  the  solution  for  velocity  Is  demon- 
strated In  the  supplement  to  this  section. 


74. 


A few  further  points  involved  In  the  failure  of  the 
method  are  discussed  In  the  next  section  with  particular  refer- 
ence to  the  problem  of  the  Incompressible  flow  past  an  elliptic 
cylinder. 

Su^/~ lament  to  section  1;  Improvement  of  the  solution  due 

~ to  coordinate  shift 

In  the  following  Is  shown  the  Improvement  in  the 
solution  for  the  velocity  near  the  leading  edge  (at  x • 0) 
effected  hy  a constant  shift  In  the  coordinate  system.  The 
airfoil  Is  assumed  to  have  the  profile 

y = €F(x)  =.  e 4-  ...]  , 0 X ^ c.  (3.8) 


If  a usual  perturbation  solution  Is  found  for  the 
problem  and  Is  expressed  In  terms  of  new  variables  X and  Y, 
the  streamline  ■ 0 becomes, 


Y4.|xffnil 


(3.9) 


or,  for  the  profile  Itself 

^ = 1 


(X-w)^+  Y^ 


(3.10) 


Now  since  the  profile  given  by  (3.10)  differs  from  that  given 
by  (3.8),  an  error  In  velocity  will  appear  so  that  at  corres- 
ponding points,  on  the  actual  profile  (3.8)  and  on  the 
solution  (3.10)  there  will  be  velocity  discrepancy  of  size 


In 


- 

(3.11) 


(velocity) 

so  long  es  the  rate  of  change  of  velocity 


11) 
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2 

Is  of  order  onej  a difference  in  profiles  of  order  £ will 

o 

introduce  a velocity  error  of  order  £~  only.  However  near 
the  ler/llni'  edge  as  was  shown  earlier  in  the  section  * 


0(— -(velocity) ) 


0( 


■2>  2-vl;(l) 

) - 0(£-— - 


2 

so  that  for  a boundary  difference  of  order  £ (with  s = 1/2) , 


0(:|^(Teloclty))  . 


0(-^)  . 


and  from  (3.11)  the  boundaries  must  agree  through  order  £ 
to  give  accuracy  even  up  to  order  £.  in  velocity. 

Thus  It  Ic  to  be  shown  that  a coordinate  correction 
of  the  type, 


X = X ^ 


+ 


= X + 6 


2 

IT 


(3.12) 


y = y 


a y 


v/lll  bring  the  actual  profile  and  the  solution  profile  Into 

2 

agreement  at  least  through  order  £ In  the  neighborhood  of 
the  leading  edge.  Away  from  this  region,  since  the  velocity 

2 

derivative  becomes  of  order  one,  a constant  shift  of  order  £ 

2 

will  give  only  an  £ -order  error  in  velocity. 

Near  the  leading  edge  on  the  actual  airfoil  let 
there  be  a parameter  k such  that 


and  then 


(3.13) 
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Then  the  same  coordinates  are  to  hold  through  order  for 


corresponding  points  on  the  streamline  of  (3.10).  In  the 
Integral  of  this  equation,  letting  Y = y = F^k  f ^ e_nd  w r t^, 


one  finds 

t^dt 

t^-2t^X+X^-^I§c^£^ 

0 ^ '0 


t^dt  _ 


r1c  (((c 

B<it  ‘ 


A 


t^tt 


0 


3 A. 


(3.14) 


where 


b = -2X 


B (b*>b^  - ‘■'O^ 

(b^“Vb^-l^c)^  - 4c 


X^  + ^ 


2 b-»-  '\l  b^-4c 


‘=1  = 


D r 


-4c 


(b-»*\)?^-4o)^  - 4c 


2 b -'Tb^-4c 



t^  - 


Then  since 


Yc 


Bdt 


t^+  t? 


® tan"^  W- 


and  since  t^,  are  order  YT  which  Is  order  In  the 


region  under  consideration, 

PTe 

Ddt  B 


t2-ht2“^iL2  -V7 


TT  H 


J3 


TT 


.1 
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and  thus 

2 I F, 


TT 


K5- 


o 

f'dt 


Jr 


= £ 


(3. 


•1 


From  (3.12),  in  the  region  defined  by  (3.13) 

C 

X - X u:-  = - -f). 


(3.16) 


2 

and  this  value  for  X together  with  y a Y * £ F^k  must  yield 

2 

a point  on  the  streamline  (3.10)  at  least  through  order  £ 
accurecy.  Prom  (3.15)  since  t-, , tg  = oVT  , an  error  of 

in  (3.16)  would,  cause  an  error  of  order  one  in  (3.10) » 
and  so  it  need  only  be  shown  that  (3.16)  as  given  allows  (3.10) 
to  hold  through  order  one.  Ey  straight  substitution  one  obtains 


-H  Ik)  - k(-f  - Ik) 


1 


as  required. 

In  other  v:ords  ne?r  the  leading  edge  the  boundaries 
differ  by  0(  £-^)  only;  away  from  this  region  this  discrepancy 
is  a.6  most  0(  £ ),  and  the  velocity  error  has  been  smoothed 
to  be  at  most  0(  £ ) In  the  entire  vicinity. 

Section  2;  Incompressible  flow  past  an  elliptic  cylinder 

In  the  previous  section  we  saw  that  each  new  per- 
turbation function  introduced  an  increase  in  order 

of  singularity  near  the  leaning  edge.  One  is  led  to  Inquire 
then  If  It  would  not  be  possible  to  use  only  the  terms 
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and  to  throw  the  solution  correction  entirely  Into  a perturba- 
tion In  both  X and  y coordinates . In  fact  It  may  be  shown 
that  the  perturbation  functions  y^^^(X,X)  always  appear  In 
the  equations  associated  with  the  respective  functions, 

BO  that  either  set  of  functions  may  be  used  alone.  Unfortunately 
this  reciprocity  also  Implies  that  the  process  will  fall  In 
these  new  terms,  as  Indeed  It  does,  but  some  facts  of  Intei'est 
arise  from  the  Investigation. 

The  elliptic  cylinder  has  been  studied  In  order  that 
intermediate  approximate  solutions  could  be  compared  with 
the  exact  solutloii.  Fox*  the  eljlose, 

^ :=  1,  £<^1,  (3,17) 

the  exact  solution  for  the  complex  potential,  F,  Is 


F = 


1 “ £ 
where  z 
F 


z - 


£ f z^  - 1 


]■ 


(3.18) 


= X + ly 


= If  + itl 


potential  function 


U = 


stream  function 
uniform  flow  at  Infinity, 


Then,  attempting  a perturbation  solution, 

F - f'  £ F^-'  (x,y) , 

an  approximate  solution  may  be  found  for  the  complex  potential 

1 


F = AU 


-*•  £(2  - )|22  _ i)j 


(3.19) 
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v^.ere  A is  r conf-tnnt  dspendAnp  on  anj'  ccr.le  cbriir'es  ent, Br- 
in.:, rncl  -..here  Z has  beer,  ufed  In  the  aenee  tL.-.t  the  capital 
letter  coordinates  were  used  In  section  1. 

The  streamline  arising  from  (3.19)  as  representing 
the  airfoil  Is 

^2  , Y^(l-I-  £ n -¥2.)^ 

A.  “r  ^ . 

^ 1^ 


Thus  when  the  (x,jO -plane  and  the  (X,Y) -plane  are 
superimposed,  a discrepancy  shows  up  between  the  airfoil 
streamline  In  the  two  planes. 

The  question  to  be  Investigated  Is  whether  a coor- 
dinate perturbation  of  the  form 

z - Z + £z^^hz)  e ...  (3.21) 


can  be  used  to  resolve  ths  difference  between  the  two  stream- 
lines and  thus  make  (3-19)  a solution  of  the  problem.  For 


the  ellipses  of  (3.17)  and  (5.20)  the  mapolng  Is  given  by 

1 


z r 


V 1+  26 


Z(l+  £-£“)+  £^-  V'Z 


(3.22) 


and  since  this  can  be  expressed  In  the  form  of  equation  (3.21) 

z sr  Z+  £^(-  I - i))  - - 1 £^(gz-t|yz^  -T)-h... 

(3.23) 

one  is  led  to  hope  that  the  perturbation  Idea  will  work  for 
the  problem. 
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In  order  to  find  the  solution  (3»23)  In  the  general 
case  one  would  take  the  perturbations, 


X = X + fcx^^^(X,y)  £^x^^^(X,Y)  -f  ... 

y = Y+  £y^^^(X,Y)  €^y^^hx,Y)  -t  . . . 

and  find  from  the  equations  that  x^^^  and  are  harmonic 

conjugate  functions  whose  boundary  conditions  are  to  be  set 
by  a perturbation  mapping  of  (3.17)  Into  (3,20).  As  far  as 
the  first  order  perturbation  functions,  x'^^  and  y'^^,  are 
concerned.  It  may  easily  be  shown  correct  to  take 


(3,2^) 


The  ecuatlon  for  the  agreement  of  boundaries  then  becomes, 

^ £^y^^^^  + 2 E^y'^^y^^^+  ...J  " 2£(i-x^)-»- 

£^(1+  3(1-X^))  + . ..  j a 

ja  £ (1-X^)  - ^li3(l-x^+  2Xx^^^-by'^^^j  4 . . . , 


(3.25) 


and  of  course  another  condition  must  bo  specified  to  determine 
both  and  For/$a?e  (1-X^)  ^ 0(  £^)  one  obtains 

y^^'  (3.26) 

,2\  f o \ ( 

but  for  the  case  (1-X  ) a 0(  £ -)  by  expanding  both  x'*"^  and 
( 2) 

y'  ' as  power  series  around  an  extremum  of  the  strcr j-.llne 
ellipse  (3.?.0)  one  obtains,  in  terms  of  the  local  variables 
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l^^e 

Tinr 

•2 


the  Gerles  for  harmonic  conjugate  functions  In  the  fonra 


(2) 


= ^00  ^ *-10^  ^20  ^ “ ^20  '*1^  ®'30  ^ • 


" ^10 '»\'^  “ ''30  V'*’  “• 

/ p \ f ^ 

Then  by  Imposing  regularity  on  x'  ''  and  y''^'  at 

Y 

infinity  (in  the  -%  scale)  one  finds  of  course  thf.t  only  a 

tr 

constant  fulfills  the  requirements,  so  that  exactly  as  in 
section  1,  the  required  coordinate  change  Is 

X^^^  a - I 

^ (3.2 

= 0. 

However,  the  disagreement  of  thin  local  condition 

(3.28)  with  the  condition  (3.26)  shows  tnat  a general  boundary 

(2) 

condition  has  not  been  found,  acid  of  course  the  functions  x' 
and  y^^^  are  not  determined  until  their  values  along  the  entire 
boundary  are  given.  That  such  a general  boundary  condition 
cpcinot  be  found  by  the  perturbation  method  we  are  using  Is 
unfortunately  true  es  can  be  shown  in  the  following  way. 

The  correct  perturbation  function  of  second  order 
is  given  by  the  known  solution  of  (3.23)  as 
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When  this  is  evaluated  on  the  boundary  given  by  (3.20)  one 
finds  for  example 


f 2^ 

and  CO  it  Is  shown  that  the  unique  correct  function  j 
which  is  specified  uniquely  by  the  boundary  condition  (3.30) 
must  have  that  boundary  condilion  a function  of  f , Any  per- 
turbation method  which  requires  the  conditions  to  be  free 
from  £ cannot  give  correct  result. 

In  suinmiLry,  then,  the  attempt  to  improve  an  approximate 
solution  near  the  leading  edge  of  the  airfoil  falls  beyond  the 
first  step  since  the  boundary  conditions  on  the  perturbation 
functions  must  Involve  £ which  contradicts  the  perturbation 
method  Itself.  Although  it  is  true  that  the  one-step  cor- 
rection involving  a constant  shift  can  reduce  the  large  error 
in  the  velocity  at  the  leading  edge,  no  clear  cut  method  for 
extending  the  process  seems  to  present  Itself. 

OiB  further  attempt  at  coordinate  perturbation 
was  attempted  and  will  be  described  briefly.  It  seemed 
poselble  that  the  difficulty  with  the  boundary  conditions 
nesr  the  leading  edge  was  due  to  the  "sharpness"  or  rapid 
change  of  slope  in  this  region,  and  could  be  avoided  by  mapping 
the  streajjillne  of  the  aporoxiisate  solution  into  a noro  nearly 
circular  form.  The  process  was  carried  out,  but  proved  un- 
successful from  the  perturbation  point  of  view. 
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The  nteps  followed  vere  flret  the  determination 
of  the  uiapplng  of  the  approximate  streamline  (3.20)  Into  a 
circle  In  the  w plane.  In  this  step  the  mapping 


Z . |(w  + J) 

(3.31) 

takes  (3.20)  Into 

a circle 

H 

, .2 

a 2 ft  :z  • t m • 

Then  the  image  of 

the  actual  airfoil  (3.1?)  under 

this  mapping 

(3.31)  was  found. 

5 Q 

For  w - ^e  ' 3 ■’?)  becomes 

a 

0 = 

'^£“+sin^Q'^  £ slnQ 

(3.32) 

r~  2 ~2. 

-it  ^ sir.  0 - e slnO 


The  question  then  la  raised  as  to  whether  ^ may  be 
expressed  as  a perturbation  series 

^ - 1 £g(Q)  + ...  (3.33) 

where  g(0)  and  Its  derivatives  are  of  order  one  near  the  region 
In  question  (0  0) . 

However 

2 £^cosO 

sln^e  aln^e  - £sln0  ) 

' (3.34) 

and  near  0*0,  0(1)  so 

2pp0  = 0(1) 

fe  = 0(1), 

and 

^=*4(^-1)  =|o(l)  = 0(i).  (3.35) 
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Thus  the  derivatives  of  g change  order  and  an  expansion  of 
the  type  (3.33)  le  not  possible  even  In  these  terms.  Aotuaily 
g Is  a function  of  t In  this  region,  and  In  fact  If  one  sets 


sin  0 


= o 


(3.36) 


then  (3.32)  may  be  expressed  as 

2 o« 


= 1 + 5> 


2(- 


OC 


l+< 


2 

=5?)  "*• 


2(- 


Cl 


3 

?) 


(3.37) 


where 
0 ^ c 


the  function 


t 


roijalns  0{  £ ) for  the  full  range 


It  seems  clerr  ■*^hat  the  method  of  coordinate  per- 
turbation as  it  has  been  used  In  this  study  la  not  satisfactory 
as  it  stands  for  improving  sc?rations  to  elHlptlo  problems. 

Clearly  If  one  takes  a coordinate  perturbation  of  the  form 

t = Z -h  £g^^^(Z)  -h  £^g^^^(Z)  (3.38) 


under  the  requirements  that 

(1)  §§  = 1 constant  as  Z"^  “» 

(11)  g*^^(Z)  be  regular  outside  the  body 

(ill)  g^^^(Z)  be  of  order  1 In  some  small  region 

(say  of  order  £ ) near  the  body  (Z  s 0(  £.  ) ) , 


then 


g'-'(Z)  r Z(cq^-~^ -h -|  + ...) 


.(1) 

.(1) 


Tr- 


i- 


or 


= function  of  -»  which  contradicts  an  exoanslon 


Z s- 

of  the  perturbation  type 
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Unfortunately  the  study  of  the  Inconpresslble  case, 
which  was  undertsJcen  In  the  hope  that  some  method  of  solution 
Inmrovenent  mlaht  be  found  which  could  ceo'ry  over  to  the 
com]  roBSlble  case,  has  not  yl<*ldecl  any  such  method.  The  in- 
corapres Bible  case  retains  its  former  status  as  a straight 
mapping  problem,  and  it  will  bo  found  that  even  the  constant 
coordinate  shift  vrhlch  served  as  a one~stage  velocity  correction 
in  this  work  cannot  be  applied  satisfactorily  to  the  com- 
pressible case. 

Section  3!  The  oompresslble  cer.e 

In  the  hope  of  finding  at  least  a one-stage  correction 
to  the  velocity  given  by  the  usual  perturbation  solution  to 
the  oompresslble  flow  past  a thin  airfoil,  the  problem  was 
investigated  again  for  the  case  of  an  elliptic  cylinder. 

For  compreeslble  flow  however  the  ea_uations  for  the  differ- 
ent order  coordinate  perturbation  functions  are  no  longer 
simple  Laplace  equations,  but  involve  nonhomogeneous  terms 
based  on  lower  order  functions.  Further,  these  terms  con- 
tribute increasing  orders  of  singularity  near  the  leading  edge 
of  the  airfoil.  Thus  where  in  the  incompressible  case  the 
coordinate  perturbation  functions  were  chosen  to  suppress 
singularities  arising  from  the  boundary  conditions,  in  the 
compressible  case  the  choice  is  governed  by  suppression  of 
singularities  in  the  equations.  In  this  case  it  will  be 
found  Just  as  before  that  the  increase  in  order  of  singularity 
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of  the  perturbation  functicne  themeelveB  prevents  the  per- 
turbation aeries  from  converging  near  the  leading  edge,  but 
even  worse  In  this  case  not  even  a first  stage  oonst*int  shift 
is  permissible  Bluce  It  does  not  satisfy  the  equations  Involved. 

That  this  la  so  may  be  seen  by  a study  of  the  equations 
of  the  problem.  The  equation  for  the  streamfunotlon,  for 
compressible  flow  la 


V I 


\ 

«2..2  i 

Mu  j 


(3.39) 


where 

- ?v  = 

M ■ Mach  number 

Cq  9 speed  of  sound  at  infinity 
^ s density. 

The  problem  studied  here  has  been  simplified  to  the  extent 
of  using  the  special  case 


2.  2 


pressure,  P = o< 
so  that  (3.39)  becomes 


- V"  = conet-nnta,  (3.^) 


yy  i 


I. 


- $ X-)  . 


^ / 
Then  solvlnr  the  problem  on  the  ellipse  (3.1?) 
2 . 

X + X -t 
£2  - 


(3.^1) 


by  the  usual  perturbation  method,  and  then  expressing  the  solu- 
tion in  terms  of  new  variables,  X and  Y,  exactly  as  was  done 
for  the  incompressible  case,  one  finds  that  a solution 
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(3.^2) 


gives  the  following  streamline  as  an  approximation  to  the 
ellipse, 


£2 


^ ^ (1  + £k)"^  r: 


where 


14- 


£ 2y2 

14  2tk 


(3.^3) 


To  correct  the  solution  a coordinate  perturbation  of  the  type 


X s:  X 


... 


Is  attempted.  (Only  one  of  the  functions,  ^ ‘ ^ or  y^^^  need 
be  used  because  they  are  redundant.  We  have  taken  ^ 0, 
1 > 1.) 


Substitution  of  (3.^4)  Into  (3.^1)  yields  a set  of 
equations  based  on  equating  to  zero  the  coefficients  of  each 
power  of  £,  . The  equations  are 


Terms  of  order  Z 

0 which  is  satisfied  by  the 
usual  first  order  pertur- 
bation function 


(3.45) 


(3.46) 
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XX 


- 3(X  - - ^,y<2) 


+x 


(1) 

YY 


-♦•  X 


(1) 

XY 


oM?.  / >b(l)^^. 
" ^ ' X ^ 

\ 


(3.47) 


Squatlona  (3.^5)»  (3.^)»  «ml  (3.^7)  ere  to  bo  used 
to  Investigate  the  poeelbls  ImprovcTisnt  In  the  solution  by  a 
constant  shift  In  the  coordinate  syBlein;  that  Is  x'^'  » constant. 

In  these  equations  the  order  of  the  derivatives  of  near 

the  leading  edge  Is  found  from  the  expression 


(1) 


- kY  - 


- X^  + A^  1-^ 


1-2(X^-Y^k^)  4-  (X^+k^Y^)  ^ 

(3.J^8) 


which  Is  the  solution  to  {3.^*-5)- 


Ner.r  Y-**0  one  flndc: 

Y(i)  ^ 

j. 


0 (3.49) 

y(  ‘o  1 

(x2-l)3/2 

t(l)  ^ 

YY  (x--l)3/2 


(?) 

Now  since  the  correction  x ' ■ conetrjit  Is  under 

consideration,  and  since  the  correction  of  order  £ , 

Is  to  he  zero  on  the  boundary,  then  as  before  we  take  x^^^  H 0 

(2' 

and  (3.^1^)  becomes  an  equation  for-  y ', 

(r  2Mir>'l)MXl)_  ^L,■(1)  Vvl)  1 .QX 

YY  - ^^2  [ ^ X ' XY  • Y ~ XX  j * (3. 50) 


This  Is  the  usual  eauntlon  obtained  In  the  ordinary 
perturbation  method  for  the  second  order  correction.  Near  the 
leading  edge  equations  (3.49)  ehow  that  the  nonhomogenous  term 
has  the  behavior 

^ or.  for  X^-l  =^e^ 

^ . (3.51) 

Neglecting  the  solution  for  y^^^  for  the  moment,  the  substitution 

x^^^  5 constant  (3*52) 


Into  equation  (3.4?)  gives 
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?' 

_ ^(2)  >u(l)  , 

J yy  ^ Y ' 


(l)x2 


)-  - (y 


+ ■'t'x '^'y]  ■ <3-53> 


In  this  equation  the  term  near  the  leading  edge 

*1  T 

ha.8  the  behavior  — y--?/?  ^ and  thus  equation  (3.^7)  becomee 

(X 

of  equal  importance  with  equation  (3.^5)  this  region  since 


3 1 

£ times  order  — r 


£ times  order 


order 


1 

£2- 

(3.5^) 


/ o ^ 

In  other  words  the  possibility  of  setting  x equal  to  a 
constant  has  been  denied  since  it  leads  to  requiring  a solution 
of  equation  (3.53),  and  It  Is  not  possible  to  effect  a simple 


correction  In  the  compressible  esse. 


♦ In  this  connection  we  mention  that  although  the  second  order 
corrections  to  the  velocities  In  the  usual  perturbp.tlon  solu- 
tions are  no  more  singular  near  the  leading  edge  then  the 
first  order  terms  for  the  case  of  Incoirpree si oie  flovr  (c.t\ 
p.  201  Clll)»  foi*  tro  c ''Tig  res  sib  le  cane  a'  Increase  In 
singularity  does  occur  at  the  second  order  level  (p.  325  C8] ) . 
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